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dome Flow Problems of Newtonian And Non-Newtonian Fluids 

The present thesis consists of the following two parts* 

Part I * Flow problems of non-Newtonian (power-law and 
Visco-elastic) fluids 

Part II * Flow and heat-transfer problem of Newtonian fluid 
in curved channels. 

In the first part (chapters II to VI ) we are concerned 
with the motion of three types of non -Newtonian fluids viz 

(i) power-law fluids, 

(ii) special type of non-Newtonian fluid exhibiting both 
power-law and Newtonian behaviour, 

(ill) Oldroyd’s model of visco-elastic fluid. 

In the second part (chapters ¥11 to IX), we discuss the 
stream line flow of a Newtonian fluid through a curved 
annulus and heat-transfer in both a curved pipe and a 
curved annulus. 

A short account of the developments of the research 
on general fluid flow problems of non-Newtoni an fluids and 
of the flow problems of Newtonian fluids in curved channels 
has been given in the introductory first chapter of this 
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thesis* Apart fro® the general survey of the allied problems* 
this chapter also contains the discussion of these flow 
problems that have led us to formulate our own research 
problems* fhis chapter also contains the discussion on the 
constitutive equations of tiiose fluids which we have used 
in our research problems. 

In the second chapter we have developed the theory 
of * ’Similar solutions of the boundary layer equations for 
power- law fluids* * for both the pseudoplastic and dilatant 
models* We obtain in this manner the generalisations of 
the Pilker-Sfean equation, Ihe velocity distribution of the 
potential flow is found to be proportional to a power of 
the length of arc measured along the wall from the stagnation 
point* Important particular cases that we have considered 
here include the following i 

(1) the boundary-layer flows along a wedge, 

(ii) the boundary-layer flows along a Hat plate, 

(iii) the boundary-layer flows in a convergent channel 

(iv) two-dimensional stagnation point flows* 

In the third chapter we have discussed the pressure drop 
in the inlet length of a circular pipe for the flow of a 
non-Newtonian power-law fluid. Following the Earman-Pohl hausen 
method, by assuming a cubic and two fourth degree velocity 
profiles, all the constants occuring in the velocity profile 
equations can be determined from the boundary conditions 
except the one * *a* * • ro determine this constant Bogus (19o9), 


XX 


on kinetic energy considerations, had suggested this * "a* * 
to be equal to h±i where n is the flow behaviour index 
of the power* law model. We have suggested and used an alter* 
native principle viz. a combination of the principle of the 
least squares and the principle of equal integrated kinetic 
energy to evaluate this constant, auad have compared our 
results with those of Bogus (!9o9)* 

In the fourth chapter we have used a sixth degree 
velocity profile using Karsnan-Pohlhausen method for the 
study of the &xially»sy metric flow of an 0 ldroyd model 
of single* relaxation-time visco-elastic liquid near a stag- 
nation point. It is found that both the boundary-layer 
thickness and the skin friction decrease with the increase 
of non-dimensional visco-elastic parameter A • the two- 
dimensional stagnation point flow has also been considered 
for the same model of visco-elastic fluid and here we find 
that the boundary layer thickness steadily increases as the 
visco-elastic parameter ^ increases. The dissipation 
function for the same fluid has also been found to agree with 
the earlier standard result. 

In the fifth chapter, the steady flows of a visco- 
elastic fluid exhibiting only the relaxation time phenomena, 
in an annulus and between two parallel plates have been 
considered* It is assumed that the fluid is injected at one 
of the boundaries at the same rate as it is withdrawn at the 
other and one of the plates (or cylinders) is moving parallel 
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to itself with a constant velocity. We have also fouad the 
solutions for the case when the velocity of injection and 
suction has any value but the relaxation time parameter is 
small md vice-versa. 

In the sixth chapter, we have obtained the exact 
analytical solutions for the flow of a non-Newtonian fluid 
characterised by the rheological equation $ 

I I . 3 ^ hd s 

t tj = Vr s 0 - 4-^a \ (_ z: z= a i 

both between two parallel plates and in a circular pipe 
under a constant pressure gradient. Ih* solutions for the 
general value of w. being very complicated algebraically, 

4<nr 

we have also given the series solutions for the fluids^ which 
either the Newtonian or the power-law term may be regarded 
as a perturbation over the other more dominant term. we have 
also obtained the solutions for some particular values of w • 
la the seventh chapter, we have carried out theoretical 
calculations for the steady, laminar, incompressible flow of 
an ordinary viscous Newtonian fluid in a curved annulus. 
Corresponding to stream function , pressure gradient \= 
and the velocity uo- for the case of a straight annulus we 
have found following Sean (=' and co' respectively 

for the curved annulus on the supposition of secondary flows 
developed in the flow region due to the curvature which has 
been assumed to be small. At first the flow parameters are 
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evaluated for the general c ase and then a particular example 
is solved numerically. The graphical representations of the 
flow line in the plane of symmetry and the projection of 
stream lines on a normal c ross-section are given. Only the 
case of large radius of curvature of the annulus is conside- 
red} more precisely this mean that the ratio of the radius 
of the outer curved pipe to that of the circle in which the 
common axis of the two curved pipes is coiled is sufficiently 
small* 

In tiie last two chapters (chapter VIII and IX), we 
shall he pre entlng our theoretical studies in the field of 
convective heat transfer problem in curved channels. These 
curved channels are * <i) curved pipe with circular cross- 
section (ii) curved annular space between two circular 
cylinder. 

In the eighth chapter we have considered the convective 
heat-transfer flow of incompressible Newtonian liquid under 
a constant pressure gradient through a circular pipe whose 
axis forms a circle of large radius it in comparison with 
the radius of the cross-section* The temperature of the wall 
is kept constant both with respect to time as well as with 
respect to direction, and both the velocity and the temperature 
profiles are fully developed. We have adopted the method 
of Imposing a perturbation temperature over the known results 
for the corresponding heat-transfer problem for a straight 
circular pipe* A number of graphs have bean drawn to show the 
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temperature profile inside the pipe and to show the local 

heat-transfer rate on the wall of the pipe giving the effect 

of the curvature of the pipe* Vie have also calculated the 

skin friction on both the convex and the concave side of 

the pipe and found out that the higher temperature difference 

in the liquid on both sides of the vertical diameter is more 

due to the value of the Prandtl number Px than due to the 

rise in skin friction, because the skin friction is higher 

on the convex side of the pipe than on the concave side* 

In the ninth chapter we have considered a problem 

similar to that of the eighth chapter i*e* here the curved 

charnel is the annular space between two concentric circles 

rotated romd an axis in the same plane distant 1 fro® 

the centre* We have considered the fully developed (both 

hydrodynasticaLly and thermally) laminar forced flow* the 

geometry of the annulus is the sane as chat in chapter ¥11, 

, u 

and the method of approach the problem A the same as that 
in the preceding chapter. Both the walls of the annulus 
(for the sake of simplicity) are kept at the same temperature* 
We have been able to evaluate she temperature profile through 
the annular region and also the local heat transfer rate at 
the inner and the outer walls of the annulus* 

the work presented in this thesis forms part of the 
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la'£t& lUwflON 

1*1 1 a4 iktfL* *Xiv il # 

fue present cmpter aim to gi ve brief account of 
the problems of Non-Newton i an fluid flows and of Newtonian 
fluid flows in curved channels, so as to place our own 
contributions to these problems in their proper perspective. 

•Section 1. 2 discusses briefly what non-Newtonian 
fluids are while section 1.3 gives some of the constitutive 
equations for these fluids, specially for those fluids with 
which we are going to be concerned in part I of this thesis. 
Section 1.4 gives an account of the work done on boundary 
layer flows of power-law fluids acid of our own contributions 
to the same which are presented in detail in chapters 11, 111 
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and IV of this thesis* Section 1*6 gives a corresponding account 
for visco-elastic fluid flows which we have discussed in 
chapter V of this thesis. Section l.$ gives a brief account 
of our investigations on non-Kewtonian fluids which are not 
included in the main thesis. Sections 1.7 arid 1.8 give a 
survey of the experimental and theoretical studies an fluid 
flow arid heat transfer in curved channels and give the stain 
results obtained by us. 

1.2 NOM-ttEtffOtfXAM FLuIDd 

fi»e rheological equation of state (or the constitutive 
equation) of a body is defined in general as the relation 
between kinematic and dynamical variables for the body. For 
simple one-dimensional flow of the simplest type of fluid 
commonly known as Newtonian fluid this is defined by a direct 
proportionality between shear-stress %. and the induced rate 
of shear ( 3u / ^ ) so that 

T = H- - ^ (1*9.1) 

** 

the const aut of proportionality, ^ is called the 
fluid viscosity. While many of the commonly known fluids ar© ,■ 
Newtonian, yet titer© are many others which deviate from this 
model. All the fluids that donut conform to the simple 
law (1.2.1) are classified as non-Newtonian} obviously, from 
mathematical point of view, the definition of non-Newtonian 
behaviour includes a far wider class of rheologl cal relations 
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than that of Newtonian fluid. Cor, sonly known non -K e wto rd an 
fluids are the slurries* soap solutions* poster* tars* various 
suspensions, hi gh polymers, greases* mineral oils* paints 
and. varnishes and a large number of lubricating oils etc* 

Non-Newtonian fluids may be classified roughly 
into the following three classes accord in - as the non- 
Newtonian viscosity y\ depends only on the state of shear 
or it depends on the previous history of the flow as well as 
on the state of shear or that the fluid exibits both the 
viscous and the elastic feature s. 

1. Purely Viscous fluids 

2 m time-dependent fluids 

3. Visco-elastic fluids. 

shaal in this thesis* discuss the flow behaviour 
of she first and the third types of the non-Newtonian fluids 
only. 

PUhilV VldwOdo fl Ul x)& 

Purely viscous fluids art those for which. the 
shearing stress is a function of rate of strain only* this 
class of fluids may further be subdivided into t 

(a) Bingham plastic Fluids 

(b) Pseudo-plastic Fluids 

(c) Dilatant Fluids 

depending on the nature of relationship of shear-stress and the 
rate of shear. Bingham plastic fluid is the simplest type of 



non-Newtonian fluids in she sense that the relationship between 
the she^r stress and the shear-rate differs from that of a 
Newtonian fluid only by the fact that, for one-dimensional 
flow, though the relationship is linear, its graph does not 
pass through the origin, thus a finite shearing stress 
is necessary to initiate the movement* the plastic flow of 
an isotropic fluid is described by Qldroyd (I9b0) by tile 
following syst em of equati ns i 


i - it 1 -*) 

* ~ \ -V- 

T = -)s§ \V -V W 


( 1 * 1 . 9 ) 


where t\ k are the de viator ic stress tensor exponents, t\ 0 
is trie const m% reciprocal mobility and cr the yield value of 
the material, equation a (1.9*9) are valid only when Jj > y <r * 
When Jj o~ * the material behaves as an elastic solid. 
Examples of the fluids which have been stated to have 
approximately Bin plastic behaviour are drilling muds, 
suspensions of chalk, grains, rocks, toothpastes, oils and 





6 


paints and sewage sludges etc* 

Pseudoplastic fluids are second in importance only 
to Newtonian, fluids* For one-d im «n s ion al flows they display 
the concave downward flow curve relationship (fig. \-l). For 
the general case the rheological equation for a pseudoplastic 
power-law fluid is given by the equation t 

1 . ' / 3 3 

W ( r. r (i.3.3) 

11 \ P.A -> I / 

wnere K and vx are called the consistency and flow 
behaviour indices respectively aid V\ < 1 * k pseudoplastic 

fluid may be thought of as a fluid compost'd of solid and 
liquid particles, which does not dilate at rapid shear rates* 
fhe asymmetric particles of its solid content, at nigher 
rates of shear align themselves along their major axes in the 
direction of the flow causing thereby a reduction of shear- 
stress as the shear-rate increases. By increasing the shear 
rate ultimately a situation may arise when no further 
alignment along the streamlines is possible. ■ fhen the fluid 
will cease to behave as pseudoplastic. It be naves thereafter 
as Newtonian fluid. 

liquation (1.2.3) gives a dilatant power-law fluid 
when vx > i • iftlatancy occurs most frequently in 
suspensions of solids at high solid content. It signifies the 
quality of dilatation as was observed by Osborne Bey no Ids who 
first perceived the chemical justification why the stress in 



7 

such fluids increases fast with the she r-rate. Ihes® fluids, 
in their state of rest, have the least voidage between their 
solid particles, which remain filled with the liquid contained 
in the fluid* So when the fluids are sheared at low rates, 
the liquid lubricates the motion of one particle past the 
other and consequently, the stresses are small* But at 
higher rates of shear, the dense pricking of the solid parti- 
cles breaks and the material dilates sll *htly causing thereby 
greater voidage between the solid particles which because of 
not being fully filled by the liquid provides lesser 
lubrication to solid particles and thus greater stresses are 
experienced at the greater rates of shear* 

¥ 1 UUU — i<L A ui’l v FL ul Pw> 

4 visco-elastic fluid possesses both the elastic 
and the viscous properties* tfncn a visco-elastic fluid is 
flowing, certain amount of energy is stored up in the material 
as strain energy, while so«e energy is lost due to viscous 
dissipation. Oldroyd*# fluid, Maxwell* s fluid, flivlln- 
krickson fluids etc. are some of foe visco-elastic fluid* 
representing dilute suspensions of solid particles in viscous 
liquids or the emulsion aid suspensions of one Hewtonian 
liquid into another* Viscoelastic effects are observed in 
polymer solutions such as poly me thyl-methcrylate solutions 
and air in touline. All the visco-elastic fluids exibit the 
normal stress effects, and that these normal stress effects 
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do not necessarily depend solely on fluid elasticity, tout 
this generally is the case, arid the usual as well as the most 
outstanding and important examples of systems exhibiting 
normal stress effects are clearly elastic in nature. 

fae pioneering studies of Weissenberg demonstrated 
the ‘climbing* of elastic fluids on to a shaft rotating within 
the fluid, in opposition to tne centrifugal forces, this 
Weissenberg effect is also manifested by an axial tension in 
such a fluid; if it is sufficiently viscous to be * ‘cut* ' 

(as in the case of molten polymer), the liquid that has 
climbed up tne shaft may be observed to contract 
longitudinally* 

nm~Dzeb\jms, Fluids 

ELme-dependent non-Newtonian fluids exhibit a 
reversible change in the apparent viscosity of more complex 
fluids at constant temperature with the duration of shear* 
this definition excludes irreversible changes due to 
permanent alternations of particles or molecules within the 
fluid, and is limited to steady- state effects in contrast to 
the unsteady-state time effects associated with visco- 
elasticity* these fluids may further be subdivided into two 
classes » . 

(i ) thixotropic Fluids 

<ii) Hheopectic Fluids 

according as the shear stress decreases or Increases with time 
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when the fluid is sheared at a constant rate. However, ve 
shall not discuss the flow of this class of fluid. 

In this thesis we have mainly confined ourselves 
to the study of only those types of non-Newtonian fluids vi*. 

(i) 0 stwald-de-Waal e power-law fluid 

(ii) Special class of non-Newtonian fluid 

(iii) Oldroyd* s model of visco-elastic fluid 

We shall now discuss the constitutive equations 
for these fluids s 

i. 3 co ns n ru cm 4 rioa a 


(1 } Oatwald-de-waale Power-Law Fluid 

4 nearly endless variety of equations have been 
proposed for portrayal of purely viscous non-Hewtonian flow 
behaviour, particularly for pseudo-plastic fluids. One of the 
models of sufficient interest is the Osfcwald-de-Waale * power- 
law* model that obeys the stress, strain-rate relationship 
characterised by the rheological equation * 


» 3 3 

= K (| 

i \ v*~-\ y 


(1.3.1) 


where K and v\ are called the consistency and the flow 
behaviour indices respectively. If h < 1 , the apparent 

viscosity decreases with increasing shear stress, and the 
fluid is called the pseudoplastic power-law fluid, whereas, 
if K >1 » the apparent viscosity increases with increasing 

shear stress and the fluid is called the dilatant power-law fluid. 
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la practical cases* the relation (1*3.1), though 
eiuperical , appears to represent a wide variety of non-Newtonian 
fluids bet car than most other proposed equations and certainly 
better than any other available two-constant equation. Fbr 
this reason the power-law rheological model for-ns the major 
portion of this thesis. Obviously, the power-law model 
includes the special case when w is unity, and K Is equal 
to the Newtonian coefficient of viscosity • denerally, 

we may say that Vv is an indent: to the degree of departure 
froft Newtonian behaviour in the sense that farther vv i 8 
removed from unity, above or below, the more pronounced becomes 
the non-Newtonian characteristic. 

In process industries the dllatant fluids are much 
less common than pseudoplastic fluids but when the power-law 
is applicable, the mathematical treatment of both types may 
be more or less the sa&e, 

(ii) Special fype of Non-Newtonian Fluid; 

It is well known that no real fluid can be said to 
satisfy one particular rheological equation completely for all 
the shear rates, fhere are some liquids which are nearly 
Newtonian under one particular set of conditions but may behave 
as power-law under another set of conditions, bills fluid 
and byring-Powell fluids are some examples of such fluids* 
in chapter VI of this thesis we have considered another example 
of such fluid which obeys the constitutive equation t 
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3 a h=I . 

t tJ - =. h^s -+ h . | ( 2 : z. a ^ ij ( 1 . 3 , 3 ) 

where € f j arc the strain-rate components, are the 

a tress components; j-*- 5 and w are the constants depend- 
ing on the flow behaviour of the fluid# This model at me 
ex.tre.ne represents a Newtonian fluid when y\ equals unity 
and the coefficient of viscosity is (^4 , while at the 

other extreme it represent s a power-law fluid for - o . 
Further for the values of v\ less than unity and small rates 
of shear, we see that the second term in ( 1.3.2} dominates 
over the first whereas for large shear rates the first term 
dominates over the second, similarly for the values of w. 
greater than Unity arid small rates of shear the first term 
dominates whereas for large shear rates the second term 
dominates* fhu® the model (1*3.3) can represent a wider range 
of behaviour than the power-law model itself. 


(ill) Oldroyd's Model of Visco-elastic Fluid 


Oldroyd, dtrawbridge and foms (1961) have shorn 
that the behaviour of dilute (Ca 3 pare* nt) solution of highly 
polymerised me thyl-methae ryl a t e in an organic liquid in a 
simple shearing motion can be approximately, at sufficiently 

gmall rates of shear, represented by the stress-strain law 
proposed by Jeffereys (1939) 






(1.3.3} 


-V- 
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wh ere is the deviator ic stress- tensor, A and 

are the relaxation and retardation fciaie parameter respectively 
of the liquid and y- is the coefficient of viscosity as in 
the case of Newtonian flow, is the rate of defor- 

mation tensor arid the dot over it denotes the change and 
is the rate of stress and is equal to the expression* 


d_ 

at 




t,. ->> k 




^3 ^ 


+ \- 


$ 


k-,* 


where 



(1*3.4) 


For incompressible liquids ^k,k -o * In this generali- 
sation, m cannot avoid non-linear equations of state if we 
remove all restrictions on the rate of strain* In (1*3.3) we 
have, therefore, considered only a particular case of this 
fluid where only the relaxation time parameter A is taken 
into account* 

Gldroyd (l95o) has considered the problem of generali- 
zing ( i*3. 3) subject to the considerations * 

(i) the equations should be valid for all shear rates 

(ii) The generalized equation should be in a form which 
does not change when coordinate system is changed, so 
that we do not need to select any particular coordinate 
system In special relationship with the flow pattern. 
The translation as well as the rotation of t&e moving 


(iii) 
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element with the motion must he taxen Into account* 
He has thus found that the equation of state for 
another Isotropic Incompressible Oldroyd liquid in Cartesian 
notations can he written as 


+ >1 [^- + * ^ ^ 


t Oxv; 


where 


= 2.^ ^ £,- k ■+ A z (2' e Ji <1,3 * 6) 

Ai >\ a ft , >?i and Vi are ‘time constats * and 


Is the static or 2ero shear rate viscosity* The 
material derivatives are given by the expressions 


-Hv) 

= 2Wi 

4- 


3|>u 

Tf~ 

at 




J)€u 

= a^ij 

-t 


2 fjA 


~ 3 t 


3 


- 60 X\ t>*j + Kv 

- -e Kj -+ 


where 




3lK 

"3xj 


a^j, \ 

3>v ) 


is the vorticity tensor. 


^ <k is the Krone cktr tensor and - t Cv 

fhis form of the equation of state implies equal pressures 
normal to the flow direction in simple shear* For determining 
any flow, this has to be solved in conduction with the usual 
continuity and equilibrium equations* 
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(iv) Rivlin- Erlcksen Fluid 

Rivlin (1948a, 1949b), Rivlin and Ericksea (1965) 
have shown that if we assume that in a visco-elastic fluid 
which is isotropic in its state of rest, the stress compo- 
nents are expressible as polynomials in the gradients of 

velocity, second acceleration Ov-O acceleration at 

the point considered, the stress matrix T = [I t* • \\ may 
be expressed as a matrix polynomial in n kinematic 

matrices A 1? A * A n whose coefficients are expressible as 

polynomials in traces of products formed out of these 
matrices, the Kinematic matrices are defined by equations. 


where 


k 


(A+t) 


K re = \\ \\ (n.^1, 2- -b) (1.3.6) 

^ >i - a. -ev j 

^ k \ j , 1 +* AhA i t5 l u, , j +■ ^ J> V 

(1.3.7) 

(.X, ha - 1,1, h ) 

U\i = > - - h) 


later on Rivlin (1966) considered certain simple 
types of steady state laminar flows vis rectilinear laminar 
flow, torsional flow between two-parallel plane discs, helical 
flow in the annular region between two coaxial cylinders, for 



th* viscoelastic fluid, tor A 


h — 0 vb®a ^- > a nd 


T = 


k J + % A, -H ^ A 


* 2_ 

i- •+- Ai A x 


-* *<s(A, A* + A* Am) 4- °<c ( A? Ai_ -t- aJ A, J 


-+ ^ 7 C A, At + A*. Ai") + J 


{ 1 # 3#8 Jt 


wfaera 1 i, the _ t , 

■*"“ q-g u, D „i T _ 

»o»lal, in the ten scalar invariants fe A, , ^ A ,» 

^ ^ ^ *>' - *>K , fc.A,'A, , k A,Cand 

ta- A% * 

It i. Obvious, that Beiner-Bivlin fluid ohar.ct- 
•rlsed by the constitutive equation 


j 


- - 1> r 1- h «vj + fe <>, v « kj (1 . 3<9) 


where f*. and fj_ are polynortrt, in invariant, II 

defined by 


«Qd Hi 


II 


= i <i{ 


lir = 3" -£yj €j* € Kn ; 


(1*3,10) 


I* only a particular case of ( 1 . 3 . 8 } vher. «* = 

- °<s- - ^ 6 ~ ^ __ ^ 2 

S«pur and Sha.hi Oort < 1961) toy. ob tained th. m 

wjUoj, , «„ #n . seised the conditions under 

* ° f ^ nWd « * stained by . 
«»if©» pressure gradient. 




1*4 «)Ui» Miff -LAYER FUr-.S OF PO^ER-L'iW fLUIdS 

As early as 1926, the empirical relation (1*3.1) 
representing the rheological equation of a power-lav fluid 
was used by seott Blair (1938) but from an engineering point 
of view the recent studies by salt (1949), Rankin ( 1965) and 
Weltmaa (1965) etc* may be of greater interest* Although 
Reiner (1946) has objected to this kind of empirical 
relation not derived from any physical concepts, even then 
in recent times this relation (1*3*1) has enabled us to 
explain various flow problems (dee olkshit (1966), Kapur (1963c), 
Kapur *nd Sriv&stava, p.M, (1963), Reiner (1960) and Proudman 
(1966) etc.) including the boundary layer flows, tnough in 
fact exact analysis of the flow situation for non-N e w toni an 
fluids is possible in only very special situations. One of 
these refers to boundary layer flows of power-law fluids. 
the boundary layer theory for Newtonian fluids is well develo- 
ped* In fact it tries to find out the asymptotic solutions 
of Navler-Stokes equations for large Reynolds number. «fhen a 
liquid Is flowing over a plate, the influence of the viscosity 
of the liquid is felt in a thin layer (the so-called the 
boundary layer) near the solid surface. Outside this layer, 
velocity gradients are small, and the theory of the potential 
flow is assumed to be valid. A number of assumptions are 
:nade to simplify these equations to find their solutions. 

Boundary layer theories for power-law fluids have been 
developed by a number of workers ( Acrivos, peterson aid shah 
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( I960 ) i Kapur (1963b) «und (1962)). flogue (1969) has 

used the boundary layer equations to extend filler’s (192 5 ) 

method for discussing the inlet length for a circular pipe 

for a Newtonian fluid, Kapur (1963c) has also integrated 

the bomdary layer equations for the flow in a two-dimensional 

Jet of an Ixicoapressible pseudoplastic power-law fluid. He 

has drawn the velocity profile for the particular case when 

the flow behaviour index w is equal to 1/3# The boundary 

y 

layer thickness has been shown to vary as H+t where k. 

is the consistency index of the fluid. Kapur and Qupta (1063a) 
have integrated the boundary-layer equations to determine 
the entrance length and the bomdary layer %r ickness as a 
function of the flow behaviour index k in two-dimensional 
flow of a po war-law fluid in a straight channel « They (1963b) 
and the author have also found that the velocity profiles in 
the inlet length of a straight channel and circular pipe 
respectively can also be approximated by f urth degree velocity 
profiles. The laminar flow of power law fluids have also been 
discussed very widely in literature (Bizzel arid, dlattery (1962), 
Kapur (1963c), Bhukla (1966), drivastava (1966) and Williams 
and bird ( 1962) ) . 

fhe boundary layer equations for the two-dimensional 
flow of power-law fluid are given by 


a lib 

dx 


•d- 


* Uoo du *J 


L 
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and 


4- ^ -O 


where u. and are the eomponen ts of velocity along and 

perpendicular to the wall, U(x) is the free stream velocity 
hri 

and q) = , h and P are the viscosity coefficient 

and the density respectively, these equations cannot in 
general be integrated exactly owing to the complexity of the 

V\, 

term (,^ u /2y.) and of Uoo being an arbitrary function. 
Integration of this equation by Kara&n-Pohlhausen method is 
not satisfactory for the entire range of ^ as has been 
shown by bizzel and battery (1963). Abbot and Kline (1960), 
Kanohar (1963) and Morgen (1962) have given the methods of 
reduction by one variable in the partial differential equations 
envolving more than one independent variables. On the basis 
of such met nods i bchowal ter (i960), lshizawa ( 1963) , Ingoff 
(1964), Oka be (1963) and many more have shown that the above 
two equations admit of * ’Similar solutions* *. duch solutions 
for the boundary layer equations of Newtonian fluids are well 
known dchlichting (1960). dchowalter was the first to give 
* ’similar solutions’ ’ for boundary layer equations of the 
pseudoplastic power-law fluid. Later on fomlta (1961) 
published a similar treatment, though math# statically 
Scho waiter’s ( I960) treatment was valid for both the values of 
y\ less than unity and also for greater than unity, yet he 
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confined himself to pseudoplmstic ( ^ <i ) fluids* fhe 
reason was Wo-fald . Experimental studies have shorn th it 
power-lsw relation describes the behaviour of a mrsber of 
real pseudoplastic fluids reasonably well in certain flow 
situations* the similar statement cannot be made for dilataat 
fluids (n > 1 ) • In particular serious difficulties 

arise woen v% approaches 2 as has been pointed out by Aorivos, 
Peterson and dhah (1960). At a later date Eapur and fyagi 
( 1964) have obtained the similar solutions of the boundary 
layer equations for two-dimensional flows of a general Keiner- 
JRivlin fluid, fhey have shown that these solutions exist if j 

the free stream velocity is perpendicular to one-third power ' 

of the arc measured along the wall from the stagnation point. 
Nanda (1962) has anowu that in the ease of Seiner- Bivlin fluid 
with constant coefficients of viscosity and cross- viscosity, 
the only possible similarity solutions for three -dimai si onal 
boundary layer equations in Cartesian coordinates are « 


Case (i) ’J = q* -v b'z- , W -v bx 


— \ 
J 


(1.4*2) 


. . U = ’M , W 
Case (ii) 




Acrivos, Peterson arid 3hah (1966) have developed an asymptotic 
method for solving the laminar boundary layer equations for 

power-law fluid under conditions where the flow external to 
the boundary layer has a general form, by the method of 
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' ‘similar solutions* 1 * Brown and Stewsrtsen (1965) have also 
obtained 1 ‘similar solutions’ * in which the vorticity decays 
algebraically towards the outer limits of the viscous layer 
and that the similarity solutions with algebraic decay can be 
the limit solutions of the full boundary layer equations with 
exponential decay. 

In chapter IX we have developed similar solutions 
of the boundary layer equations for power- law fluids when the 
velocity distribution of the potential flow is proportional 
to a power of the length of the arc measured along the wall 
from the stagnation point* A number of important particular 
cases Including 

(1 ) the boundary layer flows along a wedge 

(ii) th« boundary layer flows alon*, a flat plate 

(ill) the boundary layer flows in a convergent channel 

and (lv) two-dimensional stagnation point flows 
have been given* thm treatment of this chapter was developed 
Independently of the treatment of dcho waiter (1964) and Xomita 
(1961). in fact, a later comparison between these results 
and those of achowalter revealed a ss&fchf- natical discrepancy 
in his treatment (as reported by dkellasd (1967)). this 
descreptney in the correct form for #cho waiter** results has 
also been given in the last section of this chapter. 

In chapter III while discussing the velocity profiles 
for the power-law fluid in the inlet length of a circular pipe, 
we have found that it can also be approximated by a cubic curve 
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and a fourth degree curve* In either case a free constant *a* 

0 

remains to te determined after satisfying the boundary conditions, 
and this constant *a* we have determined by the principle of 
equal integrated kinematic energy and also by the principle 
of least squares method* these values of ‘a* are fairly 
comparable to the approximation by bogus (1959). 

In chapter VI we have considered the rheological 
equation of state in which the stress is regarded as the linear 
combination of two fcerrasf one indicating Newtonian behaviour 
and the other Indicating the power-law behaviour* This 
represents Newtonian fluid at one extreme and the power-law 
at the other extreme, and It can also represent the entire 
spectrum of behaviour between these t%ro extremes. One 
particular case of this has been used by Ng and sal be 1 (1963) 
and Kapur and Gupta (1966) to discuss the lubrication flow 
for the slider bearing. we have given the exact analytical 
solution for the flow of such a fluid between two plates and 
in the straight circular pipe. We have also discussed the 
series solutions for the fluids for which the Newtonian or 
the power-law terms can be regarded as a perturbation over 
the other more dominant term. 

1.6 VldOD-KLAdUC RU1B FLO Wb 

Viscoelastic non-Newtonian flows haw been discussed 
by .Keiner (1946, i960), Hivlla ( 1948a, 1948b), Oldroyd (1060, 
(1961), 1968) and others* Certain axially-sym metric problems 
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viz. flow through a circular tube have been considered by 
Rivlin and Green ( 19-56} , Bhatnagar and Rao (19b?), Ofeerol and 
Kapur (I960) and Rcheehter (1961). Flow in an annulus between 
two c Ircul&r cylinders has been considered by Serr in (19o9), 
Gafcfca (I96l), Fredrickson (1961), Narasimhan (1961), ftathna 
ar*d daleshwarl (1962) and Pipkin (1966)* Langlois and Rivlin 
(1963) have considered the slow and steady state flow of a 
visco-elastic fluid characterised by Rivl in-Kricksen equation 
of state (1.3*8) through non-circular tubes* Sharma (1969b) 
has applied Karman-pahl hausen method in solving the ordinary 
differential equation obtained in the study of rotation of 
a plane lamina in the visco-elastic fluid characterised by 
(1*3.4)* Langlois (1963) has extended the Newtonian flow with 
negligible inertia considered by proudman (1956) and E&beraan 
(1982) between two spheres to the case of Rivlin-Iricksen 
fluid and have shown that the inclusion of the visco-elastic 
term in the equations of motion influences the flow between 
spheres In two distinct ways s Some of these terms modify the 
distribution of the rotatary flow about the axis, whereas the 
other of them contribute the secondary flow in the aerldinal 
planes. Kapur and Ooel (1961a) have obtained the conditions 
under which the rectilinear flow of a visco-elastic filvlin- 
gr ioksen fluid can be maintained under a uniform pressure 
gradient under the sarnie conditions as obtained by Rivlin and 
Green (1966) for the Reiner-Bivlin fluid. Rajeshwarl and 
Rathna (1962) have followed the method of Rrivastava,K.C.(1968) 



in applying Karman-Pohlhausen method to obtain an approximate 
solution of the flow of a visco-elastic fluid ( 1 * 3 . 4 ) near a 
stagnation point. Datta (1964) has considered the flow of a 
particular visco-elastic fluid between two parallel porous 
plates with suction and injection. Sharma ( 1969a) has also 
considered flow of the same fluid near the stagnation point. 
Leslie and fanner ( 1961) have investigated the slow flow of a 
visco-elastic liquid past a sphere. Mohan Bao ( 1962) has also 
discussed the flow of the same liquid between two rotating 
coaxial cones having the same vertex, pipkin (1963) has shown 
that the recti liner flow of non-linear visco-elastic fluid 
through straight tubes can only occur© if the tube is of 
circular cross-section or if the rheological constants satisfy 
certain relations. Fredrickson (I9€l) and Tanner (1963) have 
discussed the helical flow of visco-elastic fluids in channels 
Williams and Bird (1969a) have discussed the steady flow of 
Oldroyd* * model of visco-elastic fluid in tubes, slits and 
narrow annulli. They (williams and Bird ( 1969b) ) have also 
studied the three constant Gldroyd* s model for visco-elastic 
fluids. The flow of such fluids has also been discussed by 
Thomas and Walters (1963) in a curved pipe under a pressure 
gradient. 3riva s tava, P.K. (1963) has discussed the propaga- 
tion of small disturbances in a seal-infinite visco-elastic 
liquid due to the slow angular motion of a disc in two cases 
vis* (i) when the angular velocity of the disc lo (0 - co 0 s O) 
(Dirac Delta function) and (il) to a) ■= to 0 
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The solution, as the relaxation time parameter * tends to 
aero, has been shown to corresponds to the Newtonian ease* 

Gberoi and Kapur (i960) gave some general considerations 
regarding axi ally-symmetric non-Newtonian flows and derived 
some new solutions for tnis problem# They have used the 
method generally used for discussing the axially symmetric 
viscous and axially symmetric hydroma gne tic flows* 

In chapter ¥ we have considered the steady motion In 
an annulus and between two parallel plates of the viscoelastic 
fluid characterised by the rheological equation (1.3*4)* fbe 
fluid Is assumed to be injected at one of the boundaries at 
the same rate as it is withdrawn at the other, and one of the 
plates or the cylinders is moving with a constant velocity 
parallel to itself, solutions have been obtained for small 
velocity of suction and injection for all relaxation time 
parameter and for small relaxation times for all velocities of 
suction and injection* 

In chapter IV we have considered the axial ly-symaetric 
and two-dimensional stagnation point flows of the visco- 
elastic fluid characterised by oldroyd*s model (1*3*4) in which 
only the relaxation time parameter A is taken into account. 

W® have been able to give a sixth degree velocity profile to 
satisfy all the boundary conditions which were not fully 
satisfied by Sharna. It has been found that both the boundary 
layer thickness and the skin-friction decrease with the increase 
of the visco-elastic parameter - the results just reverse of 
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the corresponding results of Sharaa but in agreement with the 
results of Rathna and gajeshwari (,1962) for another model. 

The reasons for this difference have also been examined* The 
two-dimensional stagnation point flow have also been 
considered* 

1.6 SOME pam&M* OM hEXKftt-BXVLlil mo fflWJH-mCStiGii FLUIDS 

In a paper we (2964) used aeiner-Siviin fluid charac- 
terised by the rheological equation (1.3,9) with coefficients 
of viscosity and cross-viscosity as constants* as a lubricant 
in an externally pressurised bearing. There we found that 
both the flow rate of the lubricant and the load capacity 
decrease with the increase in the coefficient of cross- 
viscosity if the pressure in the recess Is kept the same. If. 
however* this pressure is so adjusted as to ensure the same 
flow rate as for a hewtoaian fluid, we found that the load 
capacity actually Increases with the increase in the cross- 
viscosity coefficient, m additional effect of this presence 
of the cross-viscosity is that the surface of the liquid film 
exposed to the admospherie pressure t> 0 given by are 

no longer circular cylinders though these are still the 
surfaces of revolution. Since these surfaces are concave 
outward* these are yet another explanation that these fluids 
exhibit the Meissen berg phenomena. In another paper we (1966) 
considered a number of possible flows for Reiaer-Rtviin fluids 
when the coefficient of viscosity is constant and the 
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coefficient of cross-viscosity fc is a function of 
invariants (1,3* 10). fi.ere we concluded that though it Is 
satisfying to some extent to be able to extend the exact 
solutions for Newtonian fluids to the store general Beiner- 
fiivlin fluids (see Bhafcnagar (1964)) such an extension is 
possible only for a very special class of those fluids for 
which and fi are constants. For two-dimensional flows, 
however, ft can be allowed to vary. For Poisseulle and 
Couette flows, however, it is possible to obtain solutions 
for the general Heine r- Hi vl in fluids (see Belner (i960)). 

It is, however, easy to see that in all the following cases s 
(1) convergent and divergent flows between two non- 
parallel walls, 

(ii) two-dimensional stagnation point flows, 

(iii) axi ally-symmetric stagnation point flows 
and (iv) flow about a rotating disc 

extensions are not possible even if walls are porous. 

Vta have not included the above two investigations 
in this thesis, as recently the very existence of such fluids 
has been denoted, since the dissipation function of such 
fluids can possibly be negative (see Jtaour (19b®)), which can 
never happen in physical world. 

in another investigation we discussed the use of a 
slightly non-Newtonian liquid (Bivlin and Brick sen (1966) 
and langlois (3963) ) as lubricant In a slider bearing with- 
out any side leakage. The lower sunaer surface moves with a 



constant velocity parallel to itself while the stationary 
shoe (finite in length) above it is inclined at a small angle 
<>( supports the load. The direction of the motion of the 
runner is towards the narrower face of the bearing, ife found 
tnat the equip ressure surfaces are no longer cylindrical in 
shape but still they are surfaces of revolution concave 
towards the direction of flow, and this again illustrates 
the Weis sen berg phenomena exhibited by such fluids. Xn the 
expression for the load capacity of she bearing we saw that 
the contribution due to the elastic coefficient is always 
positive so that it increases the load capacity of the 
bearing, whereas for negative coefficient of cross- viscosity 
the contribution is still positive. the frictional resistance 
on the stationary shoe also decreases with Increase in the 
elastic coefficient, and that the difference of frictional 
resistances on the two surfaces arise only due to the Newtonian 
viscous coefficient of the fluid. 

1.7 FLOW' lii QUiiV*:.!) CH&Nm£LS 

It can undoubtedly be said that the channel flow and 
heat transfer problem are frequently encountered and play a 
very important role in science and technology. Mainly 
theoretical work has been done on three duct shape cases vis. 

(I) straight channels 

(II) curved channels 

(lii) convergent or divergent channels. 
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Comparatively* a very large number of theoretical studies have 
appeared on straight channels, fhe reason is that they are 
simplest to discuss mathematically, the curved channel 
case is very complicated from the point of view of theore- 
tical handling. Consequently, a very small number of papers, 
dealing with curved channels are available. 

In the seventh and the subsequent chapter a of this 
thesis, we are interested only in the curved channel*. Here 
a curved channel Is a duct wi«>se axis is not a straight line 
(as in the case of a straight channel) but a curved line, 
whose shape and cross-sectional ared remain m changed in the 
axial direction. 

Curved channels are encountered in several branches 
of science, specially, curved channel flow and heat-transfer 
problems are of great interest to biologists* and engineers. 
For instance, several types of curved ducts occure with in 
the human body (including the ducts through which blood flows) 
and in heat exchangers (in the heat exchange system* of 
various industrial equipments, curved channels act as heating 
or cooling coils). Construction of curved passages for the 
coolant in the heat exchange system* of various power 
producing mechanical units is nov-a-days an important problem. 
Thus, curved channels case ii of great technical importance. 

It is of particular interest from the theoretical point of 

:a.i: irjKxawi t*' :x: mpa.TX.'vmmwm 

♦One strong ^evidence of the interest of the medical scientist* 
in flow* though a curved pipe was obtained when we received 
a large number of requests, from medical institutes all over 
the world, for,, our paper on the subject published in 
Applied scientific Research. 



view also, since the field of curved etonnel flow sad heat 
transfer problem presents many difficult problems. 

the theoretical study of the subject of curved channel 
flov was started by Dean (1937). This pioneer paper deals 
with a laminar, steady, constant property and fully developed 
forced flow inside a curved pipe of simply connected circular 
cross*section whose axis forms a circular arc. The necessary 
assumption made in this paper Is that the curvature of the 
channel is so small that squares and higher powers of the 
ratio of the radius of the cross-section to that of the 
radius of curvature of the arc formed by the axis can be 
neglected. The theory thus gives us an approximation of 
the first order. It is in good qualitative agreement with 
experiment of justice (1911). However, the first experimental 
study of the subject was made by Grindley and Gibson (1908). 
The theory of Dean (1937) does not admit the dependence of 
the flow rate and of the resistance coefficient on pipe 
curvature. The next paper of Dean (1928) stows these 
dependences. The theory of this piper is based on a closer 
approximation „that of Us earlier paper of 1937. This 
extension, so far as it has been carried, represents a fourth 
approximation to the problem. Dean (193?, 1938) point out 
that the representative number of the problem ( called Dean 
number) is ol & 7 /ji * obviously De-a*s studies are valid 

only for small values of this number. Dean (1928) showed 
, that the flow rate and the resistance coefficient are 
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functions of Dean number alone. Ibis was veil confirmed by 
experimental work of White (1929), As a matter of fact f the 
aim of the experimental study of White was to test the Dean’s 
studies. White (1929) was in complete agreement with Dean 
(1927,38) and also gave values of Dean number for which 
Dean* s work is valid. For a large Dean number, we have the 
first study of Adler (1934). This contains both theoretical 
and experimental investigations, fhe theory is based on 
the supposition of existence of a thin boundary layer 
along ..he wall. Adler’s resistance coefficient formula 
asymptotically approaches the observed results in very large 
Dean number range, works of White uid Adler s-ow that the 
variation of resistance coefficient with pipe curvature is 
negligible up to certain small values of Dean number. 

Experiments of haves (1930) and squire (1964) show that in 
large Dean number case, various forces are important only in 
a thin boundary layer near the wall (see Barua (1963)). Studies 
of Haves and Squire motivated Barua (1963), who assumed that 
the fluid motion in the vicinity of the wail occur within 
a thin boundary layer adhering to the wall and the fluid 
motion outside this layer is confined to the plane parallel 
to the plane of symmetry of the channel, the resi stance 
coefficient thus obtained agrees with the observation* of 
White and Adler. In the analytical works of Adler and Barua, 
some of the boundary conditions are not satisfied at the edge 
of the boundary, the velocity profile of the inside boundary 
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layer region is not joined smoothly to that of the outside 
boundary layer region and the circumferential velocity does 
not vanish at c$> -=. — T (or ). Adler and Barua do 

not give complete or approximate momentum balance analysis 
for velocity profile and introduce discontinuity of tangsia- 
tlal velocity at the ed;e of the boundary layer and, there- 
fore, they set a finite tangential velocity in the boundary 
layer at 4> ■= -^/x • these defects in the analysis of 
Adler and Barua have been removed in a recent work of Mori 
and Nakaya.ua (1965). Also, this paper deals with a heat 
transfer problem and provides experimental observe tive for 
both flow and heat transfer problems. It should be noted 
that studies of Dean, are based on perturbation technique 
(l.e. technique of successive approximations) whereas those 
of Adler (1934), Barua (1963) and Hori and Nakayama (I960) 
make use of the technique of writing no meatus and -the energy 
integral boundary layer equations and the Harman-Pohl hausen 
technique, fhe theoretical studies referred to in the fore- 
going, deal with laminar flow case only. Most of the studies 
dealing with turbulent flow case are experimental vis* 

Jeschke (1996), Nippert (1939), Richter (1930), Detra (1963), 
White (1933), cunnings (l96o) etc. Theoretical analysis of 
I to (1969) and Mori and Nakayama (1967a) are based cm the 
boundary layer concept, and the Karman-Pohlfaausen technique. 

Jure. 

The two studies 1 ' some experimental results also. 

Following Dean's technique, some extensions of Dean 
(193?) are also available, e.g. we have Jones (1960) for 
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Heine r-Rivl in fluid with constant coefficients of viscosity 
and cross-viscosity, Clegg and Powers (1963) for Bingham 
plastic fluids, Ihofflas and Walters (1963) for an elastico- 
viscous liquid, Thomas and Walters (1964) for an elaatico- 
viscous liquid and elliptic cross-section, Pathak (1966) 
for Heiner-Hivlin fluid of constant coefficients of visco- 
sities and curved annulus and chaudhary (1964) for the t 
motion of a viscous liquid in a curved pipe of the cross- 
section of the for® Y-J^©) . Our own studies, which are 
also based on Dean* s technique are given in chapter ¥11. 

In the seventh chapter we have carried out theoretical 
calculations for the steady, laminar. Incompressible flow 
of an ordinary viscous Newtonian fluid in a curved annulus. 
Corresponding to the stream function V- , pressure gradient 
J> and the velocity for the case of a straight annulus, 

f j / 

we have found (following Dean (1937)),^ > and to resp- 
ectively for the curved am ulus on the supposition of second- 
ary flows developed in the flow region due to the curvature 
which has been assumed to be small. At first the flow 
parameter are evaluated for the general case and then a 
particular example is solved numerically. We have carried 
out our calculations on X£K 7044 Computer to find out the 
velocity field inside the annular space, and the corresponding 
superposition on the velocity field due to the curvature of 
the pipe* We find that on the concave side of the horizontal 
diameter the superposed velocity stows down the fluid particles 
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while on the convex side of the pipe it gives a little 
acceleration to then. foe graphical repress! tati on of the 
flow line in the plane of symmetry m& the projection of 
the streamlines on a normal cross-section are also given. 
However, we have confined to the case of large radius of 
curvature of the pipe. 

1.8 COKVtCnVii HMX-WJiNiiFlH IN CUHVEO CHANNELS 

The curved channel heat transfer problem is of great 
technical importance, since the coils of tubes are largely 
used, for example, in connection with heating and refrige- 
ration in order to transfer heat from one fluid to another. 
Under such conditions, the local heat transfer rate distrib- 
ution on the curved channel boundary will in general be quite 
different from that on the corresponding straight channel 
boundary, thus the heat transfer problems associated with 
curved channels are of Importance for heat exchanger tech- 
nology. Curved channel heat transfer problems are also found 
in heat engines and such various very industrial equipment 
as spiral tube heat exchangers. The removal of heat by mean# 
of a coolant flowing through curved passages in gas turbines 
and chemical and nuclear reactors may be said to be an 
important problem in the field of heat transfer studies. 

Few studies on the problem are available in the 
literature. Of these, mostly they are experimental (as a 
matter of fact the very first study is an experimental one) 
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i.e. Jeschke (1926), Merkel (1937), White (1932), Haves ( 1933 ), 
£de (1961), Seban and Me Laughlin (1963), Hogers and May hew 
(1964) etc* Tli# first theoretical study in the field of 
curved channel heat transfer is that of Mori and Wakayama 
(1966)* After this, there are two more theoretical studies 
in our knowledge, namely Mori and Wakayama (1967a) and Mori 
and Wakayama (3967b)* Ihe first two of these studies of 
Mori and Wakayama are experimental also* Each of these 
three papers deals with Dean' s curved pipe geometry* Mori 
and Wakayama (1965) give theoretical and experimental 
investigation® of heat transfer associated with a laminar, 
steady, constant property fluid and fully developed (both 
hydrodynamicaliy and thermally) forced convection under the 
uniform wall Ik at flux condition, the purpose of Mori and 
Wakayama (l»67a) is to extend their 1966 paper to include 
the case of turbulent forced convection* the purpose of 
Mori and Wakayama ( 1967b) is to give theoretical solution 
of the problems for hydrodynamicaliy fully developed but 
thermally developing forced convection under the conditions 
of constant wall temperature. In each of these studies Mori 
and Wakayama, i)ean number has been assumed to be very large 
and boundary layer concept has been used both for velocity 
and temperature fields. 

In these theoretical studies, heat generation as a 
result of dissipation of mechanical energy due to viscosity 
property has not been considered for curved pipe case. 
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An - analysis of beat generation due to dissipation for curved 
pipe is given in the eighth chapter of this thesis. In the 
last chapter of this thesis we shall be discussing that last 
of our studies i.e. the analysis of convective heat transfer 
due to the viscous dissipation in curved annulus streamlines 
in the field of curved charnel flow and heat transfer which 
were proposed to be given in the second part of this thesis* 

In the eighth chapter we have considered the convective 
heat transfer flow of an incompressible Kewtonian liquid 
under a constant pressure gradient through a circular pipe 
whose axis forms a circle of large radius JH t much larger 
than the radius of the cross-section of the pipe, fhe 
temperature of the wall is kept constant both with respect 
to time as well as with respect to direction, and both the 
velocity and the temperature profiles are fully developed* 
have adopted the method of imposing a perturation temp- 
erature over the known results for corresponding heat-transfer 
problem for a straight circular pipe. We have been able to 
find the temperature profile at all points of the cross-section 
of the pipe. Even the temperature Inside the pipe attains 
a higher value than that at the wall, and this point 
goes on shifting towards the wall as we go cm increasing 
the value of R the Reynold's number and that at no point 
on the concave side of the horizontal diameter, the temperature 
is above the wall temperature, Similar Is the effect of the 
prsndtl number Pv , but not of equal percentage than the 
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Key no M’s number R * IMi has been shorn to be due to 
the curvature of the pipe as veil. We have also calculated 
the akin friction on both the convex side and the concave 
side of the pipe and found that the higher temperature 
difference is the liquid on both sides of the vertical 
diameter is more due to the value of the Prandtl number 
than due to the rise in the skin friction, because the skin 
friction is higher on the convex side of the pipe than can 
the concave side. 

The second study on the conveeti ve heat-transfer in 
curved channels contained in the second part of this thesis 
is that of the convective heat-transfer in a curved annulus. 

It was just for the sake of simplicity we have included (in 
chapter VIII) the heat-transfer in a curved pipe in between 
the two types of studies on curved annulus the later one is 
discussed in the ninth chapter. 

Here the curved annulus is of the same geometry as that 
considered in the seventh chapter, We have considered the 
fully developed (both hydrodynami call y and thermally) laminar 
forced flow, the method of approach to the problem is the 
same as adopted in the case of heat-transfer in a curved pipe 
in the proceeding chapter. The walls of the annulus for 
the sake of simplicity are kept at the same temperatures. 

We have evaluated tKa. temperature profile throughout the 
annular space, and also the local heat transfer rate at both 
the inaer and the outer walls of the annulus. 
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CBAPfhR II 

SIMILAR dOLUflljttS OF £Uis> mWDAifft LaY hit E^UASIuMS FOB 
POWEB-L&W FLUIDS 

3.1 IlvfiO^JCIIul* 

In this chapter, we discuss the theory for ** similar 
solutions'* of the boundary layer equations for power-law 
fluids on the same lines as is usually dons (cf. Schlichting 
(I960)) for Hewtonian fluids* In this maimer we obtain a 
generalisation of the Falker-Skan equation* Important parti* 
culsr cases like the 

(I) boundary layer flows along a wedge 

(ii) boundary layer flows along a flat plate 
(Hi) boundary layer flows in a convergent channel 
and (Iv) two-dimensional stagnation point flows 
have briefly been discussed* 
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2.2 BA-iXS hqUiflOHS 

fhe boundary layer equations for the two-dimensional 
flow of a povei—law fluid are* as obtained by Kapur (1963c) 


a a M_ , \3- su 



(3.2.1) 


and (2.2.2) 

where u. and is- are the components of velocity along and 
perpendicular to the wall. is the free stream velocity, o 

is the density of the fluid and 9 = ■ hjp * 

fhe equation of continuity is Integrated by intro- 
ducing the stream function which is such that 




fhe boundary layer equations then reduce to 


( 2 . 2 . 3 ) 


'd , -h _ j-yz 2-1 ±l ~ u £*- u - a_ 9 A r \ 3U ^. i a ~ ! 

a#- "£*.23. ~ax 7 ‘ ax 23. L ' (3.3.4) 

m introduce new variables and functions s 
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j(?,D - . 

L • Utx} • 3-tx.) 

where &i is the modified Reynold's number defined by 


( 2 . 2 . 6 ) 


Ri ^ , , 

9 . 


(2.2.7) 


L and U* are the reference length and velocity 
respectively, and ^o) is a suitable scaling fmction to be 
chosen later. In terms of these variables, we get 
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ssj. = ucx) . i!£ n 

*1^' t ( 2 . 2 . 11 ) 

Substituting in (2.2.1) said remembering that for 
’ * similar solutions 1 * J- o , g. ) should be independent of ^ , 
w© get after consideration simplifications, the basic equation 


as * 
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U"\ ~ V" -v p. (i-4" 1 ) = * 


(2.2.19) 


where 
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We choose Uc*> and <^-oo so that c* and /S are 
constants. For 1 » the equation (2*2.12) reduces to the 

well known Faiker-skan equations 

6 + ^ + f 3 ^ i'-f") = «> (2.2.15) 
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In general the velocity component u. will increase 
from its aero value at the wall to the value ueo at the adga 
of the boundary layer and thus in this case au / a^. would 
be non-negative. From (2.2.5) and (2.2.11) it would then 
appear that if can be chosen to be a non-negative 
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function, life can take | ' to be non-negative. Xma %#t can 
make the as sumption that f " is non-negative function and 


integrate 
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obtained. 
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solution so 
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From (2.2.13) 

and (2.2.14) 

we have 
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Integrating the above equation for 
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Also fro® (3.3.13) and (3.2.14) we get 
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which gives on integration 
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where K is a dimensionless constant. 
Solving for Vju) and <^_oo we get 
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Fro® (2.3.13) and (3.2.14) It is seen that the result 
is Independent of any common factor of and rs » as it 
can be included in ^oo • therefore as long as * -=t= o » we 
can put <k ^ l without loss of generality* Also introducing 
a new parameter wv defined by 
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where is a constant of proportionality* 

Also from (2*2.6) and (2.3.9) 
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Integration of (2*3. IS) gives 
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where c *. and c 3 are constants of proportionality. 

3.4 PAKflCUL AK CASbS 
Case I s flow Fast a wedge 

If the angle of the wedge is -n 6 , the potential 

flow is given by 
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where o< is takas to be 1 • 

Also from (2*3.10)* (2.3.9)* (2. 2*6)* (2.2.8)* (2.2*9) and 
(2.2,16) after considerable simplifications m have 
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Integration of (8.4*8) subject to (8.8.18) would 

give and and then (8*4.6) and (8.4.7) would give 

\ 

u and * We can discuss both the accelerating flow 

( v > o ) and the decelerating flow < e V . 

Case II s XwO-di«enslonal stagnation Point Plow 

In the above analysis* m put (3 = vw - & ~ 1 * 


then we get 
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Wtt thus get this stagnation point boundary layer 
though unlike th© Newtonian case, this may not satisfy th« 
complete equations of motion. 

Case 111 i Flow Past a Flat Plat© 

In this case we put 
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than 
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0 *“ (say) 


this gives us th© ease of a flat plat© with s©ro 


incidence, and we get 



47 


'l = ^ 

L ^lh+1) J 


h-vi 


(3.4.14) 


V- r 


[ 


ih-| -, r— 
Mv\-vO . V . c j nvi 


(3.4.16) 


U - U, 


4 ( -o 


(2.4.16) 


- L 


1 -v -iH+i 

^(hv-\) . $ . c J . [4-h|~l "Jc 


X ) 


h 

h-H 


(2*4*1?) 


h~a 

n 1 *// 

f 


i - f 4-4 


- c 


(2.4. 18) 


Case IV * Flow in a Convergent Ghana#! 

This ease corresponds to #( ~o , ^ - i • fhen from 

(2.2.13) we get uc^ ■ * constant, and from (2.2.14) 
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th% imsie differential equation becomes 
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Integrating this again 
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or 



where the incomplete g, — functions are defined by 

n / ^ •X. b-l -l 
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(2.4,26) or (2.4.26) enable us to plot u./ u as a function of 

« If ^>1* the tables of incomplete beta functions (see 
Pearson (1932) arid (1934)) can be used. If o<>vi ^ i » 
numerical integration las to be used. 

the graph (2-1) gives */ u as a function of 
for pseudoplastic fluid (v>- 4: ) * a Newtonian fluid (n^l) 

and a dilatant fluid ( ^ = ■%. } • 

2.6 SCa>WALXM*6 lapOJSaUiii OF dlMILaam 3DLUflON6 

Schowal ter (1960) has independently attacked the 
problem of similarity solutions in three-dimensional case of 
the flow past a flat plate where potential velocity vector 
is not perpendicular to the leading edge, and this is a much 
more restrictive result than is obtained for Newtonian fluids 
in three-dimensional flow. 

Xhe object of the present section is to point out a 
discrepancy in the mathematical development of dchowalterii 
(1960) and to give necessary corrections. 



W© Us® the notations of schowalter* s paper* Equations 
C 14), (16) and (16) of the text are 


- h-H 


- <2 U ° r t^cr" O . c £ + ' VJ° r- V LI n 
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> 3 lx ? 


1 v.-,-—! • Lf'-ff"- 1 } + vm° [ra-i] 

a’l V\— ) V_ g,-£° L J 


(u 4 )^' 1 ax° 


1 - i^! (5F.J . ^ pf"_ w°. aJU f £ 

(U°) h ~' Sz° lo 0 ^ *'~ :l "a** ( U 0 )^- 1 




i - A {[ F*\ CSV'^^V} + kL - - .X 

V ' L u J ii x (.,o\nJ- 


h-v \ a rj 


(U°) n 


1 . x <s&"- 1 ] + u° r F6-n 

Cyj°) h "' sz L J B'x.° L -J 


J _.a0° F(5 J. - rr « U % r/' 

l v \N°') VlH 3x6 (\M°)^- 0 - 3^0 ^ 'z'^° 


hJ 




kL ( 

0 ‘ 2 '' 0 n°) v '' * ax 


Iheir correct forms should be respectively. 



> WJH 5U° r _,2 ^ . V' + ' ?U C r c ' ' 


(^ 0 ) h_1 Sx° 


(_F- FF'^ 4-1 + 


l {J<=) Vv >X° 


[Ffe-il 


VW! W U 

_ 3 .4^.3 F 'L J. . £1 F F -S w - ii Fe 

^U°)*-l (U°j > ' _ ' ? 2L “ 


NCM-0-1 ._. H + , 

(R) ' ^ H F " + ^ S "J*3 x F,/ ? + - j7^x ) (2.6.1) 


.. i 3UJ 1 r p'lpp'ii"! ^ _£>. — r p^-ll 

(ju 0 )*" 1- ' S*- e 1 J lUO) t '-l 3x r L ^ “-J 


( k U °) t '" 1 g- 


J— 

(u 0 )" - ' S-ZL* 


_L_ sjlug. 

(u°) h-l ax' 


F F 


v il- m_ 

^ u o)h-t 3 z‘e 


F6 




■v 


Kj_ 

u* (u c ) v ' 



( 2 . 5 . 2 ) 


and 


1 — p uC r<s'i&<3"-ii + yl__ .£^_ w -°.rFd -ai 

L J cvm°) *-< ax-° J 




_ _i ^_ c Fs" _ _JL — 1 ^ 3 . 6<s"_ u ° - s AlF 6 ," 

a*- e • (w 0 ) 1 '-* 2^ tw°)H-\ ax.. 


i^l 


*§• 


M-i a 


2f -i [ff*v (s^n *«*j 


+ u^Twr ta * 6 - 3) 


ffctta the modified fora of the equation (1?) of the 


above text is 



T 


o 


Cu°) 


m <9U C 


Q, W° (u c ) 


\~~h % f . yfc 


Q a (V- )0 J 


3W 


a, (U c 


.% [ v ) -w 


o c M 


v-h c 
Q i K c } 
fc a z ° 


- q 7 U°(w°) 


" 3Vm° 


- {*-) 


«3U* 


r <=Vi (w°) 




!\ \ \° 3$h': 


(3.5*4) 


the requirement of proportionality between V° and ^’ c as 
obtained fro® (2.5.4) say be stated as 


k 0 vj 


(2.5.5) 


and therefore (3*5.4) 


'~ h ,n« 

(U°) iH. 

v 3x° 


- GC { K 0 U°) 


Q j. Ko ^U°) 


(U°) 33 


Oq Ko 


(0°) djh 
~ ' 3Z ° 


q 4 (K 0 u c ) Ko 


q 7 CKoU 0 j 


•=. Q % ( KoK 


2 - V ) 

= £1 L 


\ — >1 i 


Q,o K u 


where K e is a constant t and 


i-w 

1 - Q v k 0 


\— V » 


--- Qg Kc = 


a , •=. Q „ =. 1 


i-k v 

_ Q | 0 Ko ; q 4 = Qq K “ 


We also get 


Mh VJ° _ Q, ko - -- - Qj % 

1 ' Z ‘ - 1 " a *' 

- Qq Ko S ^ £ - Qs -i^~ 

2 ^^ ^ KJr ' 

Integrating (3.6.8) for u c and y separately, 

0° = j O u -+- ) 

t= 

^ + A*°)] 


, 6 )’ ^ cU 0 

If. (2-5.6) 

■5fi Wo 

-w 

( 2 . 6 .?) 

(2.6.8) 

we get 

(2*6.9) 


(S§ II* lQ) 



where 


i>5 


- - k. A = o 3 (, 5l k. A ( (3.6.11) 

and 

\ 

T = Q s i (2.5.12) 

Mere ^ and \> are arbitrary constants and f is an arbitr- 
ary function. The constant t> corresponds to of 

our discussions in the preceeding sections. 

Integrating (2*5.12) and using (2.5.5), we get 


and 


U 


W 


, Ka 

C C *^o) 

(2.5.13) 

V 


c ^ -+ N-x®) 

(2.6.14) 

f> , o a , I>Ha 

C (_ at. 4- Az c ) 

(2.5.15) 


where 


W\ 


LKfc+v) — C K-1-) 3 

for the case when ^=o arid o(^ 0 * 

If (h-’j ) =. o 5 we get the solutions 

6 t o c + Az. 0 ) 

U° - C € 


(2.5.15) 
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B(jx°+ rv* 6 j }> 
- K 0 . c . c 


(3.8.1?) 
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i 

c 


£(!*- ■+ Ax* ) 
■€. 


h-g- 
w-v i 


(3.8.18) 


In the above equations, we can regard ^ , c and c - / as 

arbitrary constants and then (3.5.1?) and (3.5.11) would 
determine to Q io • * 

from (3.8.18), we find that the potential velocity 
is in a fixed direction and from (3.5.33) and (3.5.14) wo 
find the ** similar solutions ? * are possible when the potential 
velocity is proportional to so*ae power of the distance from 
a fixed st agnat ionstr al ght line. 

By putting K 0 =■ o and A , we get the two- 
dimensional case and the corresponding equations for 

\> (k+ \) — (h-o.) - 4 . b and e*' - -1 reduce to 


U ) 


Ua 

C C*-° 3 


(2.5.19) 


which is the same as (3.3.8) 
and > 


/> 

c 


I-IKa-V^iW 
0° ) K4 ‘ 


which again is the same as (2.3.9). 

For y=u+o - o- a ) = o equation (3.5.16) reduces to 



U 


c 


fi h at ' 5 



( 2 , 5 . 20 ) 


irtiicli again is the same as the value of uu) obtained in 
equation (2,3,13)} also we get 


3-0° > 


/ 

c 


e \> 


hz2z .*.* 

h-H 


(k-^) c i— ' ^ j^bl 
e' ^ ^*0*- 


( 3 . 5 , 21 ) 


which is the sate as (2.3.14) of our discussions, thus we 
am that 3cho waiter's treatment gives result which are the 
same as ours after the mathematical error is corrected. 



CHAFM III 


flOUNDAHT-LAYna VtLQGlH PSD KLUui K>R £HE F LOW UF A POWka- 
LAW FLUID 1U tm 1ICL El LMQIH OF A CIRCULAR PIPS 

3.1 maa juenoa 

Bogue (19S0) discussed the pressure drop In the 
inlet length of a circular pipe for the flow of the non- 
Kewtonian power-law fluids. Following Karman-pohlhausen 
method, he assumed a cubic equation for the velocity profile. 
While three of the constants occur ing in the cubic equation 
could be determined from the boundary conditions, the fourth 
constant *a* was selected to make the integrated kinetic 
energy of the fully developed profile the same as that of the 
theoretical one as the excess pressure drop is mostly consumed 
In increasing the kinetic energy, this consideration gave 

•a* - where k is the flow behaviour Index of the 



power-law fluid. 

In the present chapter we have suggested and used 
an alternative principle viz. the principle of least squares 
for determining *a* ' • Me find that the estimates given fcy 
this method agree closely with those given by the principle 
of integrated kinetic energy. 

fhe cubic velocity profile coincides with the fully 
developed velocity profile for n ^ A and n - K > but 
for other values of , it can be only a very good appro- 
ximation. fhe question naturally arises whether a fourth 
degree profile will not give a better fit. But a fourth 
degree profile means one more unknown constant. Following 
the practice in boundary layer theory, we take the additional 
boundary condition as the vanishing of the second order 
derivative at the edge of the boundary layer. We get a 
velocity profile which coincides with the fully developed 
velocity profile for - 1/3 > Va. ; but for other values 
of h » it can be only a good approximation, gvan for i 
it does not coincide with the fully developed profile. For 
this reason, we investigate an alternative fourth degree 
velocity profile in which instead of the vanishing of the 
second order derivative, we insist that the velocity profile 
is such that it does reduce to the parabolic velocity profile 
for the Hewtonian fluid. 

We thus have three velocity profiles, one of the 
third degree and two of the fourth degree, each of these 


containing an unknown constant * a* • To determine each *a* » 
we have three principles via. 

1 

„ ^ 

the least squares principle 
the principle of integrated kinetic energy. 

We investigate all these possibilities in the 
present chapter. Another fourth degree velocity profile is 
obtained by combining the principles of equal integrated 
kinetic energy and the principle of least squares. We 
determine two constants occuring in this equation fro® the 
boundary conditions and the other two can, in principle, be 
determined from the two principles. In practice, however, 
we obtain two cubic equations to determine the two constants 
•a* and 'd*, as it is not easy to obtain the solution of 
these two equations* 

3.2 VfeLOCm PROFILES 


(i) 

( 11 ) 

(ill) 


(1) Cubic Velocity profile 

let the boundary layer velocity profile In the 
inlet length of a circular pipe be given by 



^Ci-A/D). 
8/ D 


-+ 


(i-k/jO 1 (i- a/j,) 5 

b "V/j,“ +c 


( 3 . 2 . 1 ) 


where 

U- m the velocity component parallel to the x-axls 

on the section at a distance x from the origin 



a«& at a distance A from the central axis of 
the pipe. 


U * the velocity in the central plug 
^ 3) * the diameter of the pipe 
b * the boundary layer thickness and is a function 
of x only, 

the constants, a, b and c occur ing in (3,3,1) are 
determined from the boundary conditions * 


U - \J 

olA" 

A = T>-S 

- o 

aJt 

A * h-S 

&K 



these give the 

values 

of *b* and * 


■( 3 , 3 , 3 ) 


lot 


, so that the resulting velocity profile is of the form * 


£ -a.^l (3. 2.3 ) 


bit 


•57 D* 


& 3 / D 


(ii) fburth degree velocity profile No. l 

We assume the velocity profile in the form 


ik a U-N.61 + b OcMl’ - c 4. a 


u 


h/ b 


S a /5 


S a /J>* 




(3*3,4) 


the constants b v c and d are determined by the boundary 


conditions i 





u « u 

at 

K - b-s 


iy. _ o 

dK ~ 

at 

h. - b-S 

(3.2.5) 

22t( 

“ ° 

at 

h. — s>- s 



On applying these conditions (3,2.4) becomes t 

£ -- * ( 3 - 3 . 6 ) 

(iii) Fourth degree velocity profile No, 2 

We feel it advisable to assume the fourth degree 
profile which must reduce to the second degree profile for 
the Newtonian case, This requires in addition to the 
boundary conditions (3,2.2) that 

C ~ d - £(<*) (3.2.7) 

This gives, 

u _ (.7-5-0) (1-A ./D)* [*-%) (Q-1) U-hibf 

U ~ ' 3/P 3 ^ S' *3/3)3 3 SS/d'h 

£ 3 . 2 . 8 ) 



(iv) Fully developed velocity profile 

fke fully developed velocity profile for the flow 
of a power- lav fluid in a circular pipe (see Wilkinson 
(I960)} is 

Ji. „ r i - O' ^ 1 

v (*-*0 L ^ j (3.2,9) 

where V is the entrance velocity* Since the flux of 
the liquid passing through the pipe remains const mt, the 
equation of continuity at any section of the pipe is given 
by 


giving 


TTD X V 


- 77 ( VJ + itt 


/l U cl Ji 

D-% 


U. 

V 


1 

1 + K 4 S/5 + 


(3.2.10) 


(3.3.11) 


fhe values of the denominator of (3.3.11) when ^ =. X> for 
various velocity profiles are as follows * 

For the third degree velocity profile given by the equation 

(3.2.3), it is = a [ Q 3 - (3.2.13) 

For the fourth degree velocity profile given by the equation 

(3.2.6), it is ■= (3.2.33) 

for the fourth degree velocity profile given by the equation 
(3.2.8), It is = 2Q .t b - 


(3.2.14) 
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3.3 LEAST SQUARES PRINCIPLE OF FINDING **a* 1 

Since the assumed velocity profiles (3,2.3), (3.3.6) 
and (3.2.8) determined with the help of the equation* 
(3.2.13), (3.2.13) and (3.3.14) at the inlet length 
in the limit should approach as near to the fully developed 
velocity profile (3.2.9) as possible, we apply the principle 
of least squares to choose *a' so that 

j (u-ao^n (3.3.1) 

o 

is minimised. Here u. is the fully developed velocity 
given by (3.2.9) and (3.2.11), and u. t is the velocity 
given by (3*2.3) or (3.2.6) or (3.3.8). 

Using this we get 

(i) for the cubic velocity profile (3.3.3) 


3(^GOoV^\ii 1 <)^h 1 '+S-SS'h •+ ‘Si ) 
( h^+ 1 S| Y) + ^ ) 


(3*3.2) 


(ii) for the fourth degree velocity profile (3.2.6) 


2 .{ I + t5“l S C= Vf*4- 2lo S"i' *1 4- 1*7 I ) 
(S‘3^oh S + * 2 .SqSh‘ L 4 - L,8S n 4- 3,7. ) 


(3.3.3) 


(ill) for the fourth degree velocity profile (3*3,8) 


{ Q.3 0.S(=°'H*+'52>m ^>1 3 + 2-C,8S !&h l +48 t is;3h4- ) 

^2. C. 1 3*=! €= Sov-)^_^- to | 3V 4 ^’4 33^gn+ (7 6 ) 


(3.3.4) 


(X- 



The following table gives values of ' a ' for various 
values of n- for all the four eases vis. (3.3.2), (3.3.3), 
(3.3.4) and that assumed by Bogus (1933). 


Table 3.1 i Values of *a* for various values ©f a for 
all velocity profiles 


1 a 1 

.:.*;U.v*w«LBir ..vea'u :v « — r -nwiw Tees' i. _ mi 
^ ct > 

wn .f v tenvmmksr k ;,mjM«iwr 

<■ a } 

i rar. w iuh in: c::;::A : t: ' nu am . : a 

n 

for cubic 

velocity 

profile 

(3.3.3) 

for fourth 
degree 
velocity 
profile 
(3.3.4) 

for fourth 
degree 
velocity 
profile 
(3.2.8) 

CL = ' l+l - 
h 


km ; wr • w wnxi. ujm. 

0*4 

2.0920 

2.4109 

3,4434 

2.3664 

0.7 

2.4371 

2.3117 

3.4248 

2.4286 

0.3 

3.2722 

3.0747 

3.3428 

2*2400 

0 .9 

2.1229 

isms 

2.1141 

2.1111 

1.0 

2.0000 

1.7264 

2.0000 

3.0006 

1.1 

1.8970 

1.4938 

1.9044 

1*9091 

1.2 

141094 

1.4808 

1.3232 

1.8333 

1.3 

1.7343 

1.3833 

1.7434 

1.7492 

1A 

1.4639 

1.2944 

1.4938 

1*7368 

1.4 

1*4114 

1.2238 

1.6394 

1.4447 


-‘ftinuni::.Mrrj^i\i<<v::v'VV rnimm %rur ::cw*C- 


The above table shows that <c - gives a good 
approximation fer cubic velocity profile as wall as for the 
second fourth degree velocity profile t but the agreements 

breakdown for the first fourth degreo velocity profile, k 
possible reason may be that the second fourth degree velocity 


profile Is quite close to the third degree velocity profile 
as both of them are assumed to be reducible to the parabolic 
velocity profile for Newtonian case. 

4*4 IMfEChUTBD KlttEGC fcNattGY PKINClPLh TON FINDING ' a' . 

To make the two velocity profiles, via. the fully 
developed velocity profile and the limiting velocity 
profile correspond to each other, it is evident, however, 
that their integrated kinetic energies must also be the 
same. lence we choose 'a' so that 


r y> , /D 3 

\ A u - a. ^ A 


(3*4.1) 


where u. is the fully developed velocity given by (3.3.9) 
and U i is the limiting theoretical velocity given by 
either (3.3.3) or (3.3.6) or (3.2.8). 

Equation (3.4.1) gives s 

(1) for the cubic velocity profile (3,3.3) 


a ( iv,-h) C s ) 5+3 ^ 


3(3m4- 1} ^ _ ( to ' *JQ +U(> a + l $3^ ) 

= V<;+3, / 


°tQ. go 
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(3.4.3) 


(11) for the fourth degree velocity profile (3.2.6) 


3(3fc*t ) 


3 

f is; \ ( 300 ^ 31 ! 2-05-4, a -MT314, \ 

W+b ) V ' tioiio J (3.4.3) 


tioiao 



(ill) for the fourth degree velocity profile (3.2,8) 


_ _ / q, ? / lolsSci 3 — ° 4 - 3q70j^oa h-SUSxT \ 

^■(■2-h+i) C'Sh-t-S) \ / (3,4.4) 


JSach of these equations (3.4.2), (3.4.3) end (3.4.4) 
being cubic in ' q ' can be solved for ' Q ’ » but Instead 
of solving for ’ Q 1 we choose 1 c '> for convenient values 
of h and with this value of 'o' we compare the degree 
of agreement of the method of least squares or with the 

case of If $ denotes the ratio of integrated 

kinetic energies between the assumed velocity profile and 
the fully developed velocity profile in each of the three 
velocity profiles the comparison is shown in the following 
table * 

(i) for the cubic velocity profile given by (3.2.3) 



I83q 4- Ski 
o 

3t3v>*t) * 

Shi-3) 


( 3.4.6) 



Table 3*3 t Comparison with all the valuta of *a* far 
equation (3.4.7) 


■***"''* 1 wm t- \ r tr re :. i... rjM^jncru*. : \jfnrnamtv ru r : :»::v? »«r» ■r»r»L* i 

^whta ^vfaea *a f ia ^ifcan »a» la 
K r h +/ given by least given by equal 

-it squares method integrated K.B. 

«jr*eiee^*MamMiaiaK::nuiw.)iiMaefM«^.aa.:v... lwui 1 mmmm wurt :.*j£ \m*xz~VkZ..wrw: .a xtt? 


0*4 

0.9992876 

0*9944840 

1.0000 

0*7 

0.9993414 

0.9947159 

1,0000 

0*8 

0,9994934 

0.9954043 

1*0000 

0*9 

0.9998344 

0.9975247 

1.0000 

1.0 

1.0000000 

1.0000000 

1.0000 

1.1 

1,0000739 

1.0027474 

1.0000 

1.8 

1.0000445 

1.0044005 

1.0000 

1.3 

0.9999853 

1.0084589 

1.0000 

1.4 

0.9998403 

1.0112432 

1.0000 

1.4 

0.9994722 

1.0139794 

1.0000 


■ '.*. a* wa,m«aj.'*iwr «wtujic-a «: yjnewJavr ic^ka’tiii: v. 


(11) for the fourth degree velocity profile (3.3*4) 
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fable * Csupitisofl with all the Taluei of *a* fif 
aquatic® (3.4.6) 


when ^ ehsu *a* la 

^ given by least 

^_5-+L squares method 

h 

;iassir7;.jini]fs :::. iuw. 


.unn..'.k. va ./ ,tr^i w:n.|sr.|f il 'i;A.fif.riA:.M.enfA'k:iw He A'* , 

^¥hKl *«’ ii 
gives. by equal 
integrated &*I» 


d «6 

0*9979884 

1*0066686 

1.000000 

0 *? 

0*9949700 

1*0134666 

1.000000 

o»s 

0*9934404 

1*0197666 

1*000000 

0.9 

0*9977 €64 

1*0368066 

1*000000 

1*0 

0*9040601 

1.0036076 

1.000000 

1*1 

0*9804336 

1*0997838 

1.000000 

1*8 

0,9769486 

1*0466103 

1*000000 

1.8 

0*9786977 

1,0610043 

1.000000 

1*4 

0*9706160 

1*0669894 

1*000000 

1*8 

0*9476790 

1.0606064 

1.000000 

;:; uu : ;: si m w mww.mh m 

• i -.»*. i */ W—% /- an » . k \ v v i 


' a mut wriii e 


Uii) for tbs fourth degree velocity profile (3.3*8) 
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Jmtel# 3*4 i Comparison with all the values ®f *a* 
for aquation <3*4*7 ) 


T ;:; a - : wrv.7iwirirr*t; »,vu uvixsriavojnaiiM*:: :\r.. tar.*** i : m? ;. 

$ when 

r ‘ a- n±L 

- 

$whsn *a* is given 
by least squares 
method 

... :».r.#e«ri.^:.;:.ii ue.i n wmvrr va::i vnv .t ^ i ::ii*n 

i/iaiii* Kcm . ' ;v i..vr 'i«: vswr.'ET. x , .aiewv' 

^whsn *a* is 
givsn by equal 
integrated K*8< 

( a v»Y! :vAM6u:tiMiW)t osvnw 

0*6 

0*8699688 

0*8966674 

1,0000 

0*7 

0*9676344 

0.9388164 

1.0000 

0*8 

0 .9872894 

0.9863108 

1*0000 

0*0 

0*9809333 

0.9798086 

1.0000 

1*0 

1.000000 

1*060000 

1.0000 

1*1 

1*0186337 

1.0174631 

1.0000 

1*3 

1*0686397 

1*0366670 

1.0000 

1.3 

1*0393686 

1*0489907 

1.0000 

1.4 

1*0488040 

1.0877399 

1*0000 

1*8 

1,0663061 

1*0681630 

1.0000 

’liara. n 


n»a„n:..;r :a*i;str:6h6^^ . 


lie observe that almost In all eases, the agr t e nwi t 
between the three Methods is quite good* Our dissuasion v 
therefore justifies to sene extent. Segue* s eholee of 

^ 4 • “• ,a *' “• BMt 


values if we libs* 



CHAPIMt IV 


juiALLY-symiiiiu imq asamMim point 

fLOU OF A C&RMLJI VltfUD-bLlumC fLUI0 
4*1 Id TAD SUCTION 

The problem of viscous flow near a stagnation 
point had been considered by Bovarth (1934) for the plane 
case and by Press ling for the axi ally- sy ame tri c ease as 
has been reported by gehlichting (I960) arid Pai (1906)* 

For the Reiner-aivlin fluid with constant coefficients 
of viscosity and cross-viscosity, the problem was solved 
by drivastava, a. C« (1968) by using the Karnan-Pohl hausen 
method employing a fourth degree velocity profile, later 
Jain (1961) Integrated the saae basic differential equation 
by extremal point collocation method on IBi 660 computer 
at Bamberg* The results of Srlv&stava end Jain regarding 



the Increase In boundary layer thickness and skin friction 
were almost reverse of each other* Again Hajeshwarl and 
Bathna (1962) integrated the equations for a more general 
non-Sewtonian viscoelastic fluid which included the Belner* 
aivlin fluid as a particular ease* by using the sixth 
degree velocity profile and Karaan-pofal hausen method. The 
results agreed with the results of Jain* This shows that 
there was nothing wrong with the use of Karman-Pohlhausen 
method as an approximate method of integration* 

A close examination shows that the use of a 
fourth degree profile by Srivastava is not justified by the 
boundary conditions imposed* In fact* he uses the vanish* 
ing of the third derivative of a certain function at the 
edge of the boundary layer to determine a constant of 
integration occur ing in his problem but does not use the 
same condition for determining the constants in the velocity 
profile* If this condition had been used, the boundary 
layer thickness could have been determined without Integra* 
ting the basic differential equation. If all the conditions 
have to be satisfied consistently, the use of sixth degree 
velocity profile is imperative. 

Sharma (l9b9a) has investigated the axially* 
symmetric flow of a single relaxation-time viscoelastic 
fluid near a stagnation point by using a fourth degree 
velocity profile and for the same reasons mentioned above 
the use of a sixth degree profile is necessary in this ease 


also* tie hereby give the results of the use of the sixth 
degree profile in section 4.3 , tie find in this that the 
boundary layer thickness and the skin friction both decrease 
with the increase of the non-dimensional viscoelastic 
parameter* the two-dimensional stagnation point flow has 
also been considered for our node! of the fluid and the 
resul t corner out to be the same as that of Hajeshwarl and 
Hathna (1963) for a different model except for the change 
in the sign of the elastic parameter* As such the detailed 
calculations for this case have not been given* 

4.3 m of this aaumwer comm nous 

Let a visco-elastic liquid impinge perpendicularly 
on the wall at z-o and flow away radially In all direct- 
ions, so that the stagnation point is at the origin of 
coordinates. Let O ,VJ be the velocity components in the 
A- and z directions in the frlctionless flow, so that 


VJ - a^/2, f VJ — _ 2Az_ 


(4.3*1 ) 


being a constant. In the boundary layer the respective 
velocity components are given by 


u = , to = (4.2*2) 

Substituting these ia the equations of motion, Sharsaa (1969a) 
gets the equation 



= ® (4.2.3) 

V\\a. 

where 9 is the kinetic coefficient of viscosity and A 
is the relaxation time parameter and prism denotes differ- 
entiation with respect to 2 . Integration of the above 
equation gives 

aff" - t (4.2.4) 

where c Is the constant of integration. 

Boundary conditions used are 

-f(o-) =. -f\o) -0 

£(£*) = a , =0 (4.3.6) 

/If Iff 

4 { %*) =- > i>^ Q - 

where ^ is the boundary layer thickness. 

£'m first two of the conditions (4,2,6) eoaae fro* 
no-slip condition of the liquid, the next two cam® from the 
physical assumptions that u u and at ttas ‘ 

edge of che boundary layer. 

fhe last condition does not have any physical 
basis, though it determines c. = _ oA • Substituting this 
value and the first two boundary conditions in the differen- 
tial equation (4.3,4) at - 0 » we get another condition 



file last of the conditions (4*3.3) has not been used by 
dnarma (l9S9a) to determine the constants in the velocity 
profile, Subgt i tut ing 


s = •> n = Ji - 

= Jq? \ X - c> v 


z_ 

S 


(4.3*7) 


the basic differential equation and the boundary conditions 
become * 

S(nS $> 4 ' - &*'V V" -V S 3 ) - > (AVV-i* ) = O (4.3.8) 

4 i Co') - cfc Co') - o 

^ fi) -o ; (o} =. - £ (4,2,9 ) 


and 


U/ 


- A 6 fi) - 2 (U (i) ] 


By using the first five of the above conditions here, the 
boundary layer velocity profile becomes i 


4oy> = 


G?4 --a £> 3 3 S - x S 


ft 


i ~f *' 


+ - u - >1 (4.2.10) 



If we use the last condition also, we get 


S(^-<oS) = A (_S +4o^- s^S 2 ) (4.2.11) 

which will determine the boundary layer thickness in terms 
of the parameter A » and we would get the result even 
without integrating the basic differential equation, this 
is, of course, not a satisfactory position, this difficulty 
can be overcome by using the boundary conditions s 

I 

CM = Co ) — o 

dHD = S , 4(i ) =o » 4>(iA -o and - o (4. 2. IS) 

Also from the differential equation (4.2.8), we get 

11 / . 3 

<4 Co) - _ & 

and thus we would in all get seven boundary conditions which 
will be sufficient to determine all the seven constants in 
the sixth degree profile 

7. 1 <f 

+ Qi -V A 

+ + (4.3.13) 

It is obvious, that the velocity profile of degree 
less then six would not do, though we can use profiles of 
higher degree which would, however, mean more complicated 



calculation* 

We May also note that the above boundary conditions 
ar# the ones we usually use when a profile of degree six 
is used* 

4.3 mmmiuknm or jjo un dahst-i ay kr m amiss aa*- 

metm 

Using the boundary conditions (4*3* 13), (4*2*33) 
is obtained as 

v . S'^-^-lo'h a 5“ 3 3 S — 2-c.S 4 

^ - vv- - n - — a * k — n 


^-loS 

\o 



Q 

S 

A& ~ 


n 


(4.3.1) 


Integrating the differential equation (4.3*8) 
between the limits 0 and l, we have 


5. s 2 L 4> 4=' + S[V3 0 ‘ 1 

+ - o 

Substituting the values of 4 cv^) and its derivatives from 
(4*3*1)! simplifying and putting s' 1 - x > we get 

<x/ > (A4 8 - ) ' 3 ‘- — (^3J-tL(0+ V32-OOX ) sc 


-V- (\8mSoo - lyioooX') = o (4*3*3) 



■Since x-zp o | the above equation reduces to 


x 2 -t 


£ <-f X — 45-2-0 oc 4- iikoo .= d 


(4.3.4) 


for the ease A - o i.e* for non-elastie Mewtonian viscous 
liquid, this is the sane equation as obtained by Bajeshwari 
and Ratfana (1962) for the ease case. 


this equation (4.3.4) has two roots 6.499ft and 
37.56116 and a negative root. The velocity profile 
corresponding to 6.4994 closely agrees with that of Frossling 
as reported by schli dating (i960). The roots of the 
equation (4.3.3) for various values of A are given in the 
following table i 

TahLe 4.1 t Boots of ths equation (4.3.3) for various values 
of A . 
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.000 

6.499ft 

3.6494 

16.5696 

8.4447 

.001 

6.4637 

8.6434 

16.4333 

8*4101 

.003 

6.4380 

3.5363 

16.2969 

8.3763 

•004 


3.6313 

16.0363 

84063 

.006 

6.3859 

3.6071 

16.7694 

8.3377 

•008 

6.2147 

2.4989 

16.4926 

8.1688 

•010 

6*1464 

2.4789 

16.2338 

8*1016 

.030 

6.8031 

2.4087 

13.9766 

7*7687 

• 100 

3.3494 

14001 

6*1397 

6.3416 

•160 

3.3697 

14663 

6.9133 

6*8896 

*200 

3*1943 

1.7873 

6.7092 

6*1819 
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the she-sr at the wall z is given by 


t Yx 1 ^ P - A . 4=%) (4*3*5) 

• Jl 

the above table shows * 

(i ) the boundary layer thickness decreases with the increase 
in A , 

(11) the shearing stress at the wall also decreases with 
th# increase in A • 

iSharma (1959a) has found that both S and Vco) 
increase with A * the reason is not so much in the degree 
of the profile used as with the slight error with the cubic 
equation which should have been 1 

3 ( i^L, - \\<\ \) x" - (SIT2. + '3Mi A) 54 


+ - S-y'i&X') - O (4.3,§) 


If the correct equation is taken the results would have been 
the same as ours* 

4,4 j^UAfXOM OF dirt »M LlMiuJ 

the equation of stream lines is 


AT 4 0^ * constant 


(4.4,1) 


or 


£ FoO 


constant 



vh$r* 


H*) = 


S +2.C 

\i r' 


fc 


4- 


o "* •*- c 

j£p" 


M 


«ndl 



^~ ,c -* s . S-ia .c-j 

7^ * * w * J 


(4 



t%hl* 4,« |ini th* value* of R?) u again* t Z 


XftfeJL* 4*2 t Value* of R^) f« diffimt value* of ? 
and A 


*avnK*m**wm«t**fliwm*iiK xoifiiwcitu*^ \mmrh. xm: ; jo^rwiirr_t: wwr jt kx vmr^ r" 

j? X * 0 A* .01 A * *16 


0*2 

0*0844 

0*0861 

0*0301 

0*4 

0.2437 

0 * 343 ? 

0.3746 

1#0 

0.4883 

0*4960 

0*6863 

1*4 

0*9404 

0*9834 

1.0147 

2*0 

1*4336 

1*4649 

1.6304 


1*9411 

1*9836 

1*8804 


-*•*-' uk vcn riKjnwvM w.wm «.. f'i:mjt.:j*v:***oi7*iiifc:*imm"iivii^ 

4*6 II0'»lfttfUI6IOMJl. Sf4QHAfX0Ji 20X27 ROW 


4«8) 


fh* remit* of the boundary layer tMetei#** are 

apparently the rewM» of the results of fgajeatmari aid 
Batlm* (1962) but the nl«l« of the vi*ce-ola*t ie fluid 
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are also different in two eases and there is no aprlori 
reason why this should give the same results* However, 
to investigate this point further we study in this 
section the two-dimensional stagnation point flow of the 
same fluid* take 

u — ^ fly) •, ^ = - t L » (4.5.1) 

so that following the earlier method, stress tensor 
components up to the first power of ha*" are 

txx. — 2 1 *- + H- A^ L /4 + + 2.*- f 

^ = ?-*■ r j" + ] (4.5.2) 

= h- ;-*f' + 

Substituting these in the equations of motion and 
eliminating the pressure, the condition of integrability 
comes out to be t 

x s. - + jy' + A^WAf* 1 ) p =° 

a 2 - l 

Making the substitutions in terms of non-dimensional 
variables end parameters i 

JT\ — 'i- ushc^e- S = 

I ~ N V ' ~s V * 

^ ^ • 4Cn) 





fhe equation* (4. 6.3) reduce to 


S I.S4 1 - Si V £.*] -> [ 4.+' , -hU" + 4”)- (4.8.6) 

which is the same as equation (2.3) of Hajeshwari and 
E&tisia (1962) except that the sign of A has changed, the 
boundary condition* are also the same, thus in this case 
the effects of the parameter is exactly the same as 
that of -A in the discussions of Eajeshwari and Rathna 
(1962). Ha are, therefore, not surprised that the effects 
of A are in opposite direction in the t vo models. 

Again the dissipation function for our model 
comes out to be 


fA.a ' 1 






A_ 

S i 


^ o C4 4- V - 4> 4>" } 4* -V 4 4 ' j 


(4.6.6) 


where 



which is again exactly the same as the corresponding 
expression obtained by Bajeghwari and Rathna (1962) except 
for the change of the sign of a • 

to determine the reasons for this difference, let 

/' 7 \ / 

tax , and denote the contributions of the elastic 



ttrsi to the model of the fluid used by Rajeshw&ri sad 
latlma (1962) and fc^ x } and tbs correspon- 

ding contributions for ur model, then we find that 

i " ' 

I f 

w = .*1 


and 


f ,i 

it t 

= - t x ^ 


(4.6*7) 


Xhe contribution of these terms to the integraMlity 
equation of our model are 




^ [_iiz ( 4 . 6 . 8 ) 


which reduces to 




(4.6,9) 


and this is negative of the corresponding contribution of 
the model used by fiajeshvarl and ttathaa (1962). Similarly 
the contributions by the elastic tern to the dissipative 
function are easily seen to be of equal magnitude but of 
opposite sign. 

Both models taken by Bajeshwari and Kathna (1962) 
and our model (when we neglect higher powers of ) are 




approximation of the more general visco-elastic fluid models, 
and we find that they agree up to this order of approximation 
except for the change of sign, of the parameter used, the 
agreement is exact for the two-dimensional case mad is very 
good for the axially-symmetric case* 



CHAP £tiA ¥ 


QK SihhSS mauM OF A VlaCO»is2L%BflC LI MO BhfWMa £*> PLaM 
BOON MM hi Ai»i> Ui M AliaUi.Ua 41 fH a&Q MiVwUf aUCIIOA 

5.1 XdXHUJUcrioi* 

la this chapter we consider a certain type of visco- 
elastic liquid whose constitutive equation is a modification 
of the Oldroyd's equation (1.3.3) and is denoted as 

-v ^ = Sj 4 - -^ij (5*1*1) 

where t-jj is the stress tensor, the strain-rate 

tensor, t V j Is the rate of stress tensor and is equal to 
the expression 
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where are the velocity components, A is the relaxation, 

time parameter and A is a scalar coorspondlng to the 
viscosity coefficient. The flow of such a fluid in an 
annulus between coaxial circular cylinders and between two 
parallel plane boundaries has been considered by Datta (I960). 
He assumed that the fluid is injected at one boundary and 
withdrawn at the other both at the same rate. He further 
assumed one of the boundaries to be moving with a constant 
velocity parallel to the axis* He also neglected the squares 
and higher powers of the relaxation time parameter. 

We have found it possible to obtain solutions for 
all values of the relaxation time parameter A provided the 
velocity of suction and injection is small. We have also 
found the solutions for the case when the velocity of suction 
and injection lias any value but the relaxation time parameter 
is small, the second case has also been considered by 
Datta (3960) i but his equation (10) does not include the 
second order derivative. We have obtained a series solution 
for this case. 

We have derived a number of properties of the flows 
from our solutions. Further, tables have been prepared to 
Illustrate the results quantitatively. 
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PAHX OME 

FLO hi BbMbm PAKftLLEL FLAXES 
6.2 BASIC KtfUAl'lCftS MiD AftALYHC dGLOfJUu 

m choose the rectangular Cartesian coordinates 
with x-axis along the copper plate which is fixed and yaxis 
perpendicular to the plates directed into the liquid, the 
flow is assumed to be in the x-direetion due to the motion 
of the lower plate and the liquid is injected at the lower 
plate and withdrawn at the upper plate (or it »ay fee 
injected at the upper plate and withdrawn at the lower 
plate), then if we assume the flow to be steady and fully 
developed and the velocity components to be the functions 
of y alone | we have t 


U*. ^ Uty.) , ^ ^ 

the boundary conditions are 


U-- o 

at 

* 

~ O 

Vj 

at 

> 

^ -L 

= _ V 

at 

> 

-■ o 

v 

at 

* 

= ~L 


together with the equation of continuity 


<6.2.l) 


( 6 . 2 . 2 ) 



everywhere 


(0.2.3) 


^ = _V 

The stress strain-rate relations are i 


."XOt 

■+ 

A 

‘^1 

1= ° 



-V 

A 

r r^. -i 

L ' a* J 

O 

(6.2.4) 

: n 

“V 

A 

] 


au \ 

3? > 


The equations of mammtwa give ; 


iili minuting t 

( o« 3*4) we get 


n 


M- 4^ 



a 1 

o 

= a + 

from 

second and 


““ H ' 


( 0 * 3 * 6 ) 


Ml ~ ^ ) + 5 ’’ 0 ' %■ ^ - f* t (6.9.6) 

This is different from Dattas* equation no. (10) and it does 
not reduce to his equation by neglecting A* . 

From the first of the equations (6.3.0) on integra- 
ting and applying the boundary condition = -V , we get s 



where 1) is the constant of Integration* 

KLiminating fro® (0.2.6) and (a, 2*6} we get 

a second order differential equation which can at once be 
integrated since it does not contain ^ explicitly* Alter- 
natively solving the second equation of (o.2.3) we get 

- A -£ (6.2.8) 

where A is a constant. 

ye can, obtain solutions for two extreme cases vis 
when A is very small or when V is very small. this 
we proceed to do in the succeeding sections. 

6.3 mmhd soiunoa roa bhall VELOuinto of ducno., abb 

lUJSJQnJk AMD FOB All Bi&AXAflOtf flMB P ABAM e/fSB 

Assuming the series solution for ( 0.2.6) in the 

form t 

H 

a = V 

VU ~-c 

along with H K (0.3.1) 

3> = K V 

and substituting (o.2*7 ) and (6.3.1) in (6.2.6) we get the 
first few terms of the differential equation as t 
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\_ ^ Uo4 ' VU ' 4 " V +v'u?4-- ) -V ( DoA- V^u. 


+■ ^ ^ 2 > "V- £>M 4 - V !)sf. . )1 X r -^- V 4 - v'*' a M ~i- . v S ^M 7 

L ^ ^ 


-* V ^ + \/ s i* Mr _+. . 7 + $vtAp r lHs + v^ -vv’ L 2M> + .v^£l3 

ax T J I L s 3- av 


*> " 


2 e? 




+ v "^+ v 4 7^-] 3 +^ r _v 


r \ > 2^ 1 a~^ l w'i <f"Ma. y 1 * 2U <3 -7 

2 ?> +V ^ + +v ^ j* 


\j~P ^ + Vu i t-v + V^u^-v yCf Mt r -H •• ^ -V* (^ ^ts -V V -V- -v 


+ ^3)l, + v s » jt .. )] = f4r^5 +v £M J+v 1 |i f .. -j (0,3.2) 

Collecting the terass independent of V i.e. for toe case 
of no suction or injection we have 


^*■0^ _ l)o‘i> 

3 i ~ £" 


( 6 . 3 . 3 ) 


where 


V = 


£ 

~u 


Integrating this equation and applying the boundary conditi- 
ons t 


u, _ u at = A 

( 6 . 3 . 4 ) 


at 


0 


Tr - a 
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w# get 


Mo 

U 




(5.3.5) 


and 3>c = 

-"t* 

fMs result agrees with the resuit obtained by Datta (i960) 
for the case of no suction or injection. 

When the squares and higher powers of V are negle- 
cted, the coefficients of v from (5.3.2) give the linear 
differential equation 


. _P u -k , . M 

a 3 - p tr ^ 


(5*3.6) 


Integrating this equation and applying the boundary condit- 
ions * 


ve get 


and 


u, =o *t l , 1 


- O 


at 


l -- ° 


ULjL 


Pi 




D l = 



(5.3.7 ) 


(5.3.8) 
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W* observe from this equation that is always 
positive and symmetric in between the plates» and has the 
maximum value at ^ i.e. aidway between the plates. 

Again collecting the coefficients of t we once 
again obtain a linear differential equation for as * 


_ D J.-1 p*y_ _ kjj\ 7 „ upV _ x 

3^ ' yr "*■ p. & 


(6.3.9) 


Again integrating this equation under the boundary conditions 


U a 



at 



we get i 


^ o at 



(t>*3e 10) 


V = VF (4.3.11) 

and ^ __ up\ upA 

Pip’ -L 

Equation (6.3.11) shows that u r is always positive In the 
upper half region while it is negative in the lower half of 
the region. 

Similarly the differential equation for *A 3 obtained 
by collecting the coefficients of v 3 is 




(5.3.12) 
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Ihe Integral of which under the boundary condition via. 


gives 


^3 - o at 


(5.3.13) 


= o 


at 


3 - * 


fV 




(0.3.33) 


and D 3 _ 0 

r-ri 

thus u 3 is always positive if a new parameter K = -^yr is 
such that 


^ < AS (5.3.3*) 

Again the coefficients of give the linear 
differential equation for as 

DU, .M + < S - 3 * 18) 

which on integration under the usual boundary conditions 


u u — o at 


i - 1 


= O 


at 


i- 


(5.3.1?) 
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gives 


Ut, 

TT 


p\ M a 


p s J ^ 1 

+ - ~- '6C>-1K^-D 


( 6 . 3 . 18 ) 


and 3). 


u - 




u - M ~k 3 


li 


7^o 


the coefficients of v~ give the following diffe- 
rential equation for u s a s 


3Us 


6 s 


Adl- - - U 1 ^ ry ( £-1 ) ( € V- 1 £ V l + A ) 
h 7*cp 


*«, 5 


V) 


pVj— £( W£)( i-^ ) -1- ~-p - ( . i_ ' 5 2r 3 (5.3.19) 


5 h 


and this on integration under the usual boundary conditions 


Us = c at 5 - i 


=. o at 


>= - 


( 5 . 3 . 20 ) 


gives 
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Us 

U 








p\ % A 


8o^ ! 


% 

1 *u P J \ ^ 

MH) +~ it *--3-) 


(5,3.21) 


and t> 


$T =. 0 


In this way the values of u fe , m , ? ... etc, may be 
obtained to any desired extent but the algebra becomes quite 
heavy as we proceed further. 

Thus the velocity profile becomes s 


a 

u 


where 


1 + ~ t * +-^- C 

_ >C - - V " R - ( ^ It + i) + 

KiP 17 34oU^ v ^ 


V 1 ' R S _ 
7io O'* 


KH) cn-*i+o-] •+ 


rVa 

2u^X 




[i ^ 03-0 --^Vy Jtn) ‘ 6 - 3 - 23 




5.4 aOLUflOH FOH iiMALL HHLA. ATIUH fltta BUS FOR ALL VaLUbd 

of tun moan of »ug non m<d irj&g no* 

the results (6.3.8) to (5.3.21) show that 
are all Enactions of X in the fora » 

u = 1*3-} -* A <${'}'> + * H ' { ^ ) (6.4.1) 

to verify the results we choose this value of 
along with 

D = Do 4- A D| +■ ^ (5.4.2) 


and substituting those values in (5.2.6) and then the value 
of thus obtained in (5.2.6) we get 


|P (_-p + j) 6 ) 4- A (-P V 4 -V 1), } + ^(-fW + It.") X 

[ -tv +A H4 ,J r v.f ") + A 2- t-ur' + V <$" ) ] + 

apv' 1 ' CM ■+ bX^' V ) 


— -v X ^ (. 2- f -f efc 1 " 1 ) 


(5.4.3) 


Collecting the terms independent of A we get the differen- 
tial equation for ®* 



97 


vtmre R - ^ and prim® denotes differentiation with 
respect to j. * 

Integrating equation (6*4.4) and applying the boun- 
dary conditions i 

I = ° 

_ VJ 

we get 

£ 

U 

and t> 0 

-Ul-e 6 -) 


at 


3- = 


at 




i- 

1- <r* 


VjR p> 


(6.4.6) 


(6.4.6) 


the expansion of this expression for in 

powers of R eases out to be 

- 3, + liU-i) + g V->- 4 H»W) 

£u-£ - 7 ^ >u-j)(,&£n+j+i) 

( 6 . 4 . 7 ) 

the differential equation for 4 (.%) is obtained 
by collecting the coefficients of A , which again cosies 
out to be the linear differential equation 



^ ^ -v ^ 4^) 


(5.4.8) 


1 __ D 

^ hk + uv ! i ^ 


Xhe solution of this equation under the boundary conditions 


4110 =. . at 


y ^ 


- 0 at 



(5*4.9) 


is found to be 


u 


and 

£>1 


- VR -f* Lw^^.n 

4U-e R ) a L J 

a _ P 

ruv^ K 


(0.4.10) 


'i. 

dial lari y collecting the coefficients of A w get a 
linear differential equation for as * 


^ + R*(„ 


■n) 


■Re" j (6.4.11) 


and its solution under the usual boundary conditions is 
found to be 
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l € Ci-e R >) { fte\ (l-e r )C4---|) ^ 


~ R 1 
+ -e ^ 


U- €*)$<!- B_^)U ^) -** ] ] 


(5.4.12) 


and t> 


oy* ft 1 * ft <P 


[ j><r R - (i-€~ R )(i-|)'] 


5.6 Dl3CU33IGh3 OF 1H£ HfidULfd 

The aquation (5.4,1) can now be written as i 


u 

u 


a-i R s- 

l-e7 R 




L R S 

0-V*)* 


« E ci-i 6 ' 1 ’) 3 c« R - ct-^) y-4.) J 


+ (l-5 e } \ 3 (6*5,1) 


where 


A ^ A v 


We have tabulated the values of ^ 4 s ) and H"^) £ a 

the whole region between two plates and observe that if 0 

the velocity is the saae as for the elastic and inelastic 
liquids (ii) for A - - • 1 we see that the existing 



velocity profile becomes more bulging forward for the values 
of R- -s due to the superposition of the terms 4 <■%) 
and 4 (v) * For increasing R the primary velocity 

field { ky ) attains the velocity of the upper moving plate 
more rapidly than for smaller R • the contribution due 
to 4<-y) goes on increasing with the increase in R but 
for values of R > M this contribution starts decreasing* 
Similarly the contribution of 4 (v) term increases till 
R - M- but starts decreasing when R attains larger 


values than 4* 
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(!*•• ^ w,W u ) «v«r the velocity field 

flvea fcf-RvV u fer various values of 
Beynolds number R end A = Q =. • 1 


ij^Mi u .vtxi jieceeimMeeerr* ASltNWMnefif^fp^etv '. »*.;*<« 

R 


,6 


1*6 


3*0 


ttirew m. wtuWt. uMnxMii 'ouaiawh mi m 


<mW 4 et© 8 *© 

.r,:ix._r*v.v:iM^ er^tmMHfic * ..." A w»u i i jt*r- ■ nuvn .,: ;, *Muwt& it 


00 

•000000 

•060000 


•000000 

.000000 

0.1 

.001307 

.00666 

.084 760 

•099327 

*286484 

0*2 

.002178 

.009220 

•038148 

.139733 

•866800 

0*3 

.003761 

.01130 

•0*3496 

.126021 

.173383 

0*4 

.00306* 

•01206 

.043273 

.107366 

.102325 

0.6 

.0030767 

.01176 

.039322 

.083996 

.066631 

0.6 

•003867 

.01066 

•033003 

.0610ft 

.089481 

0.? 

.002418 

•00864 

.086308 

•0*0818 

.014432 

041 

.001783 

.00637 

.016981 

.023980 

.006366 

0*0 

.000070 

•00336 

.008384 

•010508 

.003166 

1*0 

.000000 

.00000 

.000000 

•000000 

.000800 


news ^urur \*niu‘.eriVU i : u;kx ui i'.,:tin t wrwn*kt«<**uj. M:e«tm , .\iiwvv'(iim*.: «^iS***:r*W wKWmtirm tv :$m .mtirv mrnteiu* l*. 


Jpuato 
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PAdI II 


PLOW IN AN ANNULUS 

6.6 BASIC SqtfAXlQ«S 

We choose the cylindrical polar coordinates (h. } &,x) 
with the 2 -axis along the axis of the cylinder. In the 
steady state of the liquid flow the velocity pattern is given 
by 


U(JI) 5 = &(e) - O , - u>iA.) (5.6.1) 

If we consider the flow at a sufficient distance from the 
end so that u , \> and co may be taken to be independent 
of z and then the stress, strain-rate equations become 

A -3™ - TZ ^ kh-^ * h 

( 6 . 6 . 2 ) 

k«- -V- A L '^ ] - 0 

tsMj> -v a f -j- ~ 3 - X 

L 3A A 

“ ^<9-~2L o 


fhe equation of continuity is 



+ 7C -.° (6.6.3) 

Ihe equation of momentum is 

f ^ - au _ S ^Ax. _> ! L. 

^ ' ~dK~ * £ (6.6.4) 

I he boundary conditions are 


\> — O-c at A- =. a. 


— ^ b at /*_ -=-fc 


(5.6.5) 


where a and b are the radii of outer and inner cylinders 
respectively. 

file intier cylinder is moving parallel to its own 
axis with velocity U such that 


u -s. o at A- - c. 

=. U at A- - b 

fhe equation of continuity on integration gives 


(6.6.7 ) 


_ b = ^ (6.6.8) 

A ~ A 

fh# equation (6. 6*4) on integration gives 





10 * 


A , Lf. 


'Jh 


A- 


u. 


I> 

V 


(5.6.9) 


wh«r« 3 is the constant of integration. 

£11 sins ting (r^ from the first and second equati o n 
of (6.6.2), we get 


I 




<3M 

2X 


+ A r 

L A 


PA^ 


-JU’ZTL ■ 

Bp~ 




4>M_ L 
‘PA 



(5.6.10) 


Substituting the value of from (5.6.9) in (5.6.10), 

we get the differential equation in the fora 


f ( Ja tt u-v^ . aAa ^. 1 2U , Apa 7 --^ 2 Ao.D^ n 

PA L A A ' "a*- ’ ^ /iA A 3, j 

(6.6.11) 

transforming this into non-dimensional fora by substitutions 


M 

Vb * 

A* =. 

Mi 

Q 

t> 

i - 

A 


q 

± . 

i 

0 = 

U 

Po |\>a R 


Vi 


( 6 . 6 . 12 ) 
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we get 

<^L r , S>' . 

a-? L -f •* j* 


.A 1 -v- 3^ ’ -i 

" " jjl * ^ C't J 


-JV 


I a^ 1 f u' 


a-t" 1 - 


l l 


D ; 


■« + f 


= P lff) 


( u.6. 13) 


This differential equation is non -linear and of the second 
order and the complete analytical solution is not possible. 
We therefor®| proceed to solve (6.6.13) by the series method 
as before. 


6.7 3MLE6 ijOLUflGW tOH drfALL dUjUAaoS TIM BUf K>a ALL 

vuLocmhs of iiucnoii mo injection 

Let us assume the solution of (6.6.13) of the form 


u ' ■= ) 


l 

-V > 


-4- > V M-' C"* ) 


(6.7.15 


together with D 1 


l 

b o 


x' t>; 



D 


I 


Collecting the terms independent of >' from (6.6.13) 
after substituting (6.7,1) for the inelastic case, we get 


^ JL 1 

a* ~ * -i 


Dc 


JR 

-f 


(6.7.3) 


Xhe solution of this equation under the boundary conditions 



10 ? 


is 


4- 


Jfd) 


at 


= at 


-5-1 


-S ^ <r _ 


<5. 


u 


' Av 0 . 

t cr R - 1) 


(6.7.3) 


(6.7.4) 


and 


vj' 

(,cr*L 1) 


Similarly collecting the coefficients of A from the 
equation (6.6.13) the differential equation for a? is 
found to be of the form 


_a dp 
"<9i‘ 


4 


d 

K -i 


(j'RlRVi) -A 


(_Cr^_l) -§ ■ 


(6.7.6) 


whose integral under the boundary condition# 


*=. O 


±- O 


at 


at 


-i - i 


-1 - cr 


(6,7.6) 


is 


and 


4. 

u' 


t>. = 


RCR+^1 r_/'^ R +/it <r *_ < r R -\ f R- /I 
a L 


r cr 




“I 


H<r R -i) L a- R _l 


(6.7.7) 
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Again collecting the coefficient a of ^ from <6.6.13), 
after substituting fro® (6.7, 1 ), the differential equation 
for 4- is found to be of the fora s 

^ ^ -V- B> 


where 


- BR (R+i) ( 


i. — cr 


-S 


(6.7.8) 


u'r 

aC^-1) ^ 


Ihe solution of which under the boundary conditions i 


4- - o at 'l ^ 1 


at 


-S s <r 


(6.7.9) 


is 


U' 


: CR-^J(cr-i) V 

8^-l) * 

r ji.M D R-M „ v / , R l 

J (_ cr _ cr ) - (o" — 1 ) ~i 4- (cr-l)-£ J 


R-2 r K K-2 
SLR (cr -1 ) ^ (cr-cr ) 


R fU 


O ~ 4) 


H 

s. 



J 


(6.7.10) 



1 



Slid 


*>v = 

& tcr^. 1 ) 


[ CR-X) c 0-5-1) <> R - cr R ' H ) 


+ C^tcr R ) ] 

5.8 Sfc«lJi3 SOLUTION R>ft SMALL VILOOIflaii OF SUOIXON AMO 
INJ LOTION BUT FOB ALL BELAXAflOi HMfi 

We assume the solution of (5,6.11) of the for® 


^ - Mo 4 \?- e U, 4 . . . . 

and ( 6 . 8 . 1 ) 

^ ~ ^ !>i 4 a?-,, 1 

Substituting these values In (5.6.18) and collecting the 
terms Independent of •»>* l.t. the case of no suction or 
injection, the differential equation for u 6 Is 

. av ±! = A 0 4- (6.8.8) 

dh. p- A 

Integrating this equation under the boundary conditions 

u„ ^.o at A - Q 

(6.8.3 ) 


U at 


b 



m get 


and 


U.o 

TT 




where -{ n X 

a_ 


-tu 

U p 

jkcr 

and <T - — 


( 6 . 8 . 4 ) 


this agrees with the result for ordinary visco-elastic 
fluid for the case of no suction and Injection. 

Collecting the coefficients of -J-q, fro» (6.6.12) 
after substitutions from ( 6 . 8 . 1 ) we again get a linear 
differential equation for u, i 


< 3 ui 

xx 


I uop W SU&A 

I 4 - V. . — * w , 


^ F'W u- 


jlucr 


( 8 . 8 . 6 ) 


Integrating this equation under the boundary conditions 


M i — o at 


- o at 


A, ■=. a 


h- 


( 6 . 8 . 6 ) 


ve get 





<* t 

Jl l 



and 


U/ 

u 


ivyQ- 




r C^- 

L oi b^ 


+ 


U«.A__ 

UUV 


_Au3^ 

- ( J. _ J. 

Autr 


qp M 

Ak -h/b -j 


Ik J 

r 

or 

H 

VJQ f 

0 *— V3 ’I* 

ah 


( 6 . 8 , 7 ) 


Similarly collecting the coefficients of ^ from (6.6.1s) 
after substitutions from (6.8.1) we get 


_e>Ux _ \_ ^ tjc^Af J_ r a-b" 1 - h?_ +-L 1 

3 A- pv A PJUaO- /hL_ *- Ji^cr A* 1 ,s,v J 


r j_ ^ Q - b7 :.l . i_ 1 

” A JUa- L ^ « v fe v h> W r J 


Oq' 2 -^'"' iu (A/t) 
iocr 

(6.8,8) 

and this can also be integrated under the same boundary 
conditions vis 


u a — — o at 


-*i. :=- <L 


at 


=. o 


/*- T=^io 



113 


6.9 UlaUUSSItiK Of THE RESULTS 

Ih© solution obtained for s&all relaxation tlsi 
parameter a' can be written as 


u 

u’ 


c -/-o 

(.O' -L ) 


RUUs.) a 




*2. ~L 

RC.R-VJ-) A P Of R-4 o dc R P 

* R k-2.) C — 1 ) 1 (o' — a - ) - C 0 " — M 3 + (cr^-i.)i | 


R -i \i- 


S C<t k - 1) 


+ 2R A ) { Co- -cr^) + C<r R ”2 L) -$ — ( <r R -l) S ^ ~j 


(6.9. 1) 


where P. - 5^1 and >’ - A -— 

V? ' * 

In the following tables we have calculated functions 
•fl-s) , and M'f-O respectively. Jfrom these ire observe 

that (i) if P. - o t the velocity is the same as for 
elastic and inelastic liquid (11) if o < P. < d_ velocity 
decreases as compared with inelastic liquid, (iii) if R- d_ 
the velocity remains unchanged (iv)if R > s_ the 
velocity increases as compared with inelastic liquid. 
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8*4 i ft* aon«di*«asi©nal velocity fie!4 %■, = 

for Inelastic liquid (i.e. a’ = = o ) 
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Sabi* &«£ t fh* superposed velocity / u' over t &e 

aoa-diaeiieiaaal velocity field l Vi/ = ^Vi/ 

for •last!* liquid A = ^ - .1 end varlou* 
values of fieynolds auaber R . 
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CHAPTIB VI 


FLOW OF A SPECIAL IXPE OF NGN-NEWTQN1 AN FLUID BETWEEN TWO 
PARALLEL PLATES AND THROUGH A OIHCULAR PIPE 

6.1 INTRODUCTION 

In this chapter we have obtained the exact analyti- 
cal solutions for the flow of the liquid characterised by the 
rheological equation (1.2.2) between two parallel plates and 
in a circular pipe. The exact analytical solutions are, how- 
ever, very complicated in the general case and therefore we 
have discussed only some special cases in greater details. We 
have also discussed the series solutions for fluids for which 
either the Newtonian or the power-law term can be regarded as 
a perturbation over the other sore dominant ter®. The two 
parts of this chapter deal with two separate cases via. 
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(i) flow of this fluid between two parallel plates 
(il) flow of the same fluid in a circular pipe. 

6.2 BASIC STATIONS 

Let us consider the fluid with constitutive equation 
(1.2,2) flowing between two parallel plates, a distance *“&* 
apart under a constant pressure gradient. Using the Cartesian 
coordinates for the equations of motion when the velocity 
components are of the form : 

U U^j , ^ — o and w = o (6.2.1) 

The components of the strain-rate are given by 

—o , =- ^ ( 6 . 2 . 2 ) 

and the stress components are given by 


L_L 

- C 'Yir 


t 


^ * 


rt£) 


V\~| 


I eu_\ . ^d- 
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The equations of motion for Cartesian coordinates with constant 
pressure gradient are s 
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fhe second equations of motion (6.2.4) gives ^ =l 0 «Mch 
shows that theVe is no pressure gradient in the y-di recti on 
and the pressure profile is independent of ^ . Here 
is the pressure gradient that exists along the length of the 
channel* and the first equation of motion gives * 

3 1 3 - 3 r u. ^ * u\^ \ ^ l 

Tx " ** \_ ^ ** ^ ' ** ^ J (6.2.6) 

which on integration gives 

* 

^ %. ^ ^ l %l ' = + A (6.2.6) 

Since is essentially negative throughout the channel 

we write (6.2.6) as 
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The boundary conditions are 
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6.3 BCACT SOLUTION 

^ , 




Let us put 

3M 

/ 

- a 

and = P 

qx 

such that 

(6.2.7) is written as 
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(6.3.1) 



Differentiating this with respect to u we get 

77 • = L ^ ^ ] (- IS* ) (6.3.2) 

which on integration gives 

Pa _ -A.m. ./_ u ^ c (6.3.3) 

1 1 ~ h-*\ 1 

Also we have (6.3.1) written as, since being zero 

h ■ 

Multiplying (6.3.4) by 
we get 

h~l j d' ^ -+- lb Pu’v^ 

In principle, this equation along with the original equation 

can be solved for u' (i.e. the velocity gradient) for any n 
The constant of Integration c being determined from the 
ellminanfc of o' . Bit this becomes too much algebraically 
complicated and as such we give here the solutions for a few 
convenient values of in viz* y\ = > 1 » ^nd 2 as this 

is the range in which generally lies for the power-law 
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(6.3.4) 


3- (-u'J and adding to (6.3.3) 
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fluids. 
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Gas® 1 t n » 1/2 

The two equations (6.3.5) and (6.3.4) reduce to the 
following for the particular value of n * 1/2 


- IP o' £ 4- 6(C- Pcl) = 

p-u + u' ( /4 1 + ) + p 1 ^ 3 " - , 
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The ^ eliminant gives 
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Applying the boundary condition that u = 0 at ^ ^ , we 

get another equation for c as t 


36c. jo. £ C ( SP -f ~ Jg P h Fi -t* ) 

+ 4- *-P Po ) =o (6.3.8) 

Again elimi na ting c. from these two equations we get a 
biquadratic equation which on solution will give the velocity 
profile. 



Case 2 : n » 1 
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For n * 1, the problem reduces to the Newtonian 
case where the coefficient of viscosity is (M-h) instead 
of • 


Case 2 

to 


and 


n * 2 

For n = 2, the equations (6.3.5) and (6.3.4) reduce 


p-u 1 + 4, p u U + 4 (C- Pu) = 
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(6.3.9) 
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The u' eliminant is given by 
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P x 3- ^ h fi + 6 Pur + sPVy + ifc py«) 


- 3C e.* ^- a - Ct ( f* 3 + to -I 2. = c (6.3.10) 

and this equation under the boundary condition (6.3.8) gives 


36 c*. 6C (f<V 6f*fc?A )- PV^\- IfeP-th) = c (6.3.11) 

Again eliminating c. between these two equations we get a 

blqua dratic equation in a which on solution gives the 
velocity profile for the particular value of n » 2. We could 
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m 

try more cases for different values of n but because of the 
algebraic complications* we give herewith an approximate but 
convenient method of solution of (6,3.4) by the method of 
series. 


6.4 SEHXE3 SOLUTION mm ft IS SMALL 
The original equation 

n- U) + tu-4i-)\= - (&*+M 

under the boundary conditions (6.2.8) becomes 

-^3- = W'!*) ^ ^ C- If ) h (6.4.1) 

Assuming the series solution of the form 

- j~ = ^ t t- ' ( 6 . 4 . 2 ) 

and substituting the value of (- £M /< 5 ^) in (6.4.1) and 
collecting the coefficients of various powers of ^ , we 

get the relations s 
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and this on integration under the boundary condition that 
u. 55 0 at ^ - -L transforms into 
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Transforming the above equation (6.4.6) into non-dimensional 
form by substitutions 
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vo get 
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Again writing £ . R - ^ » we have 
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total flux © through the charnel of width -i is given by 
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which on simplification reduces to 
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6,5 SfiHIfcS SOLUTION WHIN h IS ad ALL 



Assuming th© series solution of the form 
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and substituting this value of in (6.4.1) and 

collecting the coefficients of various powers of p*- w® get 
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and th© series for [3<ya;?.)Is given by 
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Integrating the above equation under the boundary 

condition that u = 0 at y. = ve get the following velocity 
profile 
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Transforming the above equation into non-dimensional form 
as before with the help of substitutions (6.4,7) we get 



where 
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the total flux Q c when p- is small is given 

by 



where -U is the width of the channel. 

6.6 PLOW THROUGH A CIRCULAR PIPS 

Let us consider the same fluid (1.2.2) flowing 
through a circular pipe under constant pressure gradient 
along the axis of the pipe. The components of velocity are 
given by 


— 0 o and co — (6.6.1) 

The components of the strain-rate are given hy 
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The components of the stress are than given by 
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Tbe equations of motion reduce to 
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Ihe latter equation on integration reduces to 
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Since ^^j2)h) i* essentially negative as oo decreases as a 


increases. 


This equation is essentially the same as derived in 
the case of flow between two plates considered in previous 
sections of the present chap ter $ except for the chiwige of 



flwia behaviour constants ^ *»d Fk mA Tl U. pi &m 

of ^ * «nd therefore we got the sane solution in this case 
also* 

6.7 DISCU^SIOMS 

Tba series far ^ gives fey the aquation (6,4.9) 
is convergent for snail Ri and n < 1 * Iter 1 vs 
g«t i <**F>* aoi this equation reduces to 

" YT ^ ' ~ ^ ^ {6*7 *l) 

which is the parabolic Telocity profile for Newtonian fluid* 
Iter other values of n we can regard all the taros of 
(6*4*9) as the successive approximating to the velocity 
profile further epproxlna tioas being very smell* is 
illustrative exwiples we take n ■=. and k-.i t the 
equation (6.4*9) can be written as t 
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y*) _ X. |!ld V<v - . i § t&is rtteesi 



i 

X. 


(j-f) 






13 ? 


Tb# fellevin* table illustrates the successive approximations 
to the velocity profile* given by (6.7,9) and (6.7,3). 

F3f® tha so no may conclude that further ap pr oxi m atl ea i 
aill mmaia practically unaltered even if wo toko aero 
and more toms of tho velocity profile - . 

fable 6.1 i tho successive approximations to the velocity 
profile (6.4,9) for K - • 4. . 
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CHAPflH ¥11 


ifamiui flow ihbjUqh cusveo mu ulus 

7 *1 INTHUDOCnoW 

fhe problem of the flow of fluids la a curved channel is 
of great importance, both from the biologists and engineers 
point of view, m have already given in section 1.7 a brief 
account of all the investigations on flow in curved channels. 

In the present chapter, we shall extend the results of 0e*n 
(1927) for a curved pipe to the case of a curved annulus. In 
section 7,2 we shall describe the geometrical problem, la 
section 7,3 we shall give the mathematical formulation of the 
same. In the next three sections 7.4 • 7*# i* shall obtain 
expressions for ^ and f 3 which determine the seconds* 

ry flow. In section 7*7 and 7 .8 we shall give the expression 
for the stress-components. In the next three sections 





(7*9 * 7*ll> we give a numerical example and for this example 
we obtain the projection of stream lines on the normal section 
and the flow line on the central plane. In the last section 
we shall be giving the effect of radii ratio v on the velocity 
field* 

7.2 GBQHETKX OF THE P BOIL IK 

Otir annulus is the regie® between two co-axial pipes 
of circular cross-sections. The radius of the outer pipe is 
a. and that of the inner pipe ii b, The line of centres 
of the annular cross-sections is coiled in a circle with 
centre & and radius i * This line is called the central 
line* The coordinate system used is that adopted by Dean 
(1927) in discussing the curved pipe flow problem and is 
represented in figure (7.1). 0-9 is an axis through o and 
perpendicular to the plane (called the central plane) of the 
central line of the annulus, c is the centre of the annular 
cross-section by a plane passing through ov asking an angle 
<9 with a fixed plane (e « o) through ol . oc i« perpendi- 
cular to oV and is of length i • Any point P f of the 
annular cross-section ® * constant is referred to by the 
orthogonal curvilinear coordinates h-, 4> and <$ , where /i 

is the distance and 4 is the angle made by c P with 

the line through c. parallel to . The line clement is 
then given by 

tpls) l ' + (/U*) Z + U+*K4) 1 (d®) J (7*2*1) 


?*3 HAfHiJUTlOAL KJBKULAIIOfi 

In the curved pipes, there exists a secondary flow 
(Barua (1964) and others)* fhe saae will also exist in the 
curved annulus, because the fluid particles near the central 
portion of the flow region which have a higher velocity, are 
acted upon by a larger centrifugal force than the slower 
particles near the walls* Ibis leads to the emergence of a 
secondary flow which is directed out* wards <l*e* away from 
the centre of curvature) near the maximum velocity region 
and near the Inner pipe wall and is directed inwards (i.e. 
towards the centre of curvature) near the outer cylinder wall* 
30 in the light of the above physical facts, the resultant 
flow should be three dimensional. In the case of regular 
flow (i.e* laminar flow), the flow field will be symmetrical 
about the central plane (i.e. the plane whose normal is o9 
and which passes through oc • 

m proceed to discuss the steady and fully developed 
motion so that the velocity components in the coordinate 
directions A t & J <4 may be written as vH rt, 4>) ^(>, 4 ) 
and * all independent of <3 . fhe rate-of- strain 

tensor has than the physical components ; 
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fluid * 
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tiie stress tensor. 

has the physical components « 
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Hep* U3 and Xj are the stress and rate-of-sfcraln 
tensor* ^ is an undetermined isotropic pressure ? £ vj is 
Oie Kronecker tensor and p- i s the coefficient of viscosity. 
We aasuaie that th« fluid properties are physically constant 
quantities and the fluid notion is a forced notion. 

Xn absence of the body forces* the equations of linear 
momentum are 1 
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fhe equation of continuity is 
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If the annulus were straight -/ % would vanish, and (7*3.3) 
to (7.3*7) would be satisfied by the expressions t 
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together with u - o 5 v =. o , and cr = , 

We assume, for the purpose of mathematical analysis, that 
the annulus is slightly curved, that is \ and \ are 
small quantities of the first order. Following Dean (193?) 
we write 


\J - u • V - • W = oo \ 4. (aJ cwA P = b\-v V 

* ~ (7*3*9) 

where u , \> ; and \> are referred to secondary flow and 
are taken to be snail of the order of ^ and also u , 
and uo are independent of <9 • w > and I 5 , are the velocity 
and the pressure terms for the corresponding problem for 
straight annulus* We know that i 
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K = - a-z. + b 

neglecting terms of order ^ L the equation of continuity 
becomes 

an . . m. + Jl . d _± . c 

3X y- A 3 4: - (7,3*11) 

This suggests the streaa function of the fora * 
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u ^ ”a> 5 ~ aA (7.3.13) 

Substituting (7.3.9), (7.3.10) and (7.3.3) in the equations 
of notion (7,3.4) to (7,3.6) and neglecting the small teras 
and simplifying to a convenient fora we find that (7*3*4) to 
(7.3.6) reduce to t 
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Since ^ and to *?* independent of & , (7.3.15) suggest* 
that we emu take 
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The physical picture of the problea and equations (7.3.13), 
(7.3*14) and (7.3.15) suggest that * , u and ^ should 
be of the forms 
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I f , 

vhmve ^ ^ and p are functions of "1 = 4 only. 

Substituting for 4- , c*j and j> in terms of the dimensionless 
functions ^ > w ’ and ^ i a (7,3.13) to (7.3. 1*), the 
following equations are obt aine d 
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where R is the aeynolds number given by 
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Eliminating |> between (7*3.10) and (7.3.20) we get 
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7.4 aOLUXlOM FDR 

Equation (7.3,24) can be written as s 
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Differentiating each tern of (7,4.3) we obtain 
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I'US v q$ I^ti / 


_ 2 _ 

IS? 



^ -f S lu-f * C f R> ^-Ji- &>>-) -i%< 


-v- t ^ sX 4 - (MJ 


(7.4.4) 


Her® E , F , G end -H are constants of integration. 

In order to determine these four constants ve use the follow- 
ing boundary conditions t 


m- B) = 


av' 


~ o 


at -f 


(7*4.6) 


X- 4 ) 


4 ^' 


at •••£ - <t 


(7*4.6) 


i ,, 

The boundary condition =* gives 


Sf_SF + U,G^» <«y) (7.4.7) 


The boundary condition 


*»> 

■ 3 ? 


1 1 


gives 



4- t G 4- H 
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t; ^ + 4f R> g, - l^_ (say)(7.4.8) 


Tli« boundary condition H-'(a-) ^ gives 

Sl-CT £ -+■ ( 14 cr x ius _ Scs- a }F -i- 14 


{ ~ a*' - Scr^iua- + 5 lct M (imr) 7 ] 


{ 5 a- M _A. cr £ + .-^ C- & ,W _ c<o- M .Lcr £ RA 


•+(_--^a-^— jct- C )R = L 


(*ay) (7.4,9) 


And the boundary condition ^ f gives 

4 = <r 

G(T E •+- (Sct' 2 '^ | fe cr 1 - i^cr) F + I GO-' 5 ' CS - 3^F\ 


| stl a- 1 * _ llcr^^CT -v&cr (^u(r)^j R-X i 

-\- 1 U cx M _ cr *" -*■ 4^ cr^iucr -I 'Q.c-^Iact ] R-A 


-V C^a- 8 - ±S) fc - 


(say) (7.4.10) 


Solving (7.4.7) to (7.4.10) v find that 





F _ i 

. — x 

\lScr Cl~c r 1 ) ^ (j+a-'i) ^cr _t_ (__ 

[_CT Cl-o-M -!iAu o-^ L, 4-cr 2 - \ Silver + 0-<r' 1 )] k 


2lcr 


La- _ J Sa-^Lg + ] 


(7,4.11) 


V = 


L Cl-cr' 2 ) L l , 


-+• Ci-l-3cr' L ) i-- 4- a-” 7 L u 



3 1 cr *■ | ( 1+ a-' 1 ) Am <y 4- (I — cr T ) | ” 


(7.4.12) 


= k( L ' +L -) - -fc e 

^ ^ E + -t, f -v f ) 

Mow substituting ths values of E , F , G and W fr©a (7.4.11) 
t© (7.4.14) ia (7.4.3) and (7,4.4) we can completely determine 

‘-h 1 aui ‘l-— and hence u. and V- for arMirary Value 

C/ rf 


(7.4.13) 


(7.4.14) 


of radii rati© cr 
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7.5 30I.UI1OH FOB wf^ ) 

In order to calculate uj 1 we rewrite (7.3.21) in the 
following fora t 


M { -k Jq ^ Gt +4 


(7.5.1) 


substituting the value of ip 1 froa (7.4.3) in (7.6.1) an! 
then integrating we get 


cd (-<, ) - _X 


4* +Cil-|- £u +i RA<s j 


+ C-q - i * ft -fk E + ,4rr ^ )"S 


32. 


lol4 


> 1 


+ t «iii + ™ - nt RE ) * -■ ^ «'* 


Civk 8 ! aj^+ ^ rVa U-t 


- ( 4r ■ ^ + ~ R^) ^ 


3 o7%. 


- t i R p + i R'a 1 ' ♦ i rM 3 )^ 


+. ( Ij -^»5 vi RAri fL-i 


+ ^ Mu 4--^i + ^ -> (f^U 


-L- 

»>8 


(7.6.2) 



1 i<? 


vtere -T and J are constants of integration* 
fo evaluate these constants w use the following boundary 
conditions > 

/ 

co C i ) - 0 ' 

° > ^ ( <r ) - c 

the boundary condition eJ(i) r c gives 


(7*6*3) 


T -iX 


C- i F - 


+ U-lu 


S + K ^ -,5 WE + 4 




3 c? ^ 


- ( R v +-Z- r\a ) 

v IQL 9x1 C / 


-+• 


1 1 S'i 


1 . R l - 


M 


(say) (7*6*4) 


and the boundary condition fa-) -~o gives 


- j _ jla 1 = 


( 3 4--S. RT- J- RA E-4- 1 I2.V) cT^ 

V 9 k ^ Its, 31 ' 
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( -! . D r , 7 D i , S 

V \ 1Jrx K ■+• k A ) cr 


~ k^A) a- -v _ L fi* * 


^ <AU _ ( Jj ^ J1 rV 1 ' ‘•Jt, 


lun. ^ y ( 


Ci RR 


6 VJ* 


+ A RAUVcr- + X. R^' x ^ C (X juT ) 


R^A 3 cr V (jLr/ 


RA Fcr 1 (XcrJ 


(i> - ^ -tr RA <5 + -1 R A f ) cr^jL 


(say) ( 7 . 5 * 5 ) 


Solving ( 7 * 6 * 4 ) and ( 7 . 5 . 5 ) va obtain 


J* .-. ( Hr — A, J 
O-^J 


( 7 . 6 . 6 ) 




O-'crM (7,5,?) 

Thus calculating 60 we have found the axial component 
of the secondary flow, 

7.6 SOLUTION POH f>* 

(1.3,30) can be rewritten as t 

Differentiating both sides of (7,4.2) with respect to 'S 
and then substituting the expression for ^ ^ <q ^ ) in 
(7.6.1) we find 

f>’ = - -£ + Ci E +£.U + 

— Jl £-5 " — RA ^ 4- 

_ J- RA^ ^ / + U ( 1 (7-6>8) 

7.7 DiPHLdiilON fOR f© n_ 

The frictional force per unit area of the external 

boundary in the axial direction is given by 



1 Ar. a 


i 


o! 


-$72. 


* h [■ 


dtM I 

~JT 


. JL + JL(i^). 

Q + S ' > 


oo | 4 


*■] 


! 


" ^ “X‘(^ ) f _; W< ^] (7*7*1) 

and the frictional force per unit area of the inner boundary 
is given by 


^-<9 I 

1 h Hb> 



5_ / <=W 

JLA d-i 


j W<$ 
'i-tsr 


(7,7.2) 


Here “J.J' is to be calculated from (7.6.2). 

7.8 EXPHfcocilON FOR { 

The frictional force per unit area of the outer 
surface in the circumferential direction is given by 


1 ^hzXX. 

_ a^ 1 - r 7 _ j_ ^fci +. ..±1 1 css 4 

"iiT - ^ 4-s -4 1 - 

(7.8.1) 

and the frictional force per unit area of the inner surface 
in the circumferential direction is given by 



K>4, 


fat- f d^ ! 



-i ^cr 


{7 *8#2) 



Here ^ and are given by (7.4.3) and (7.4.4) reap- 

actively, also the value of can be calculated by 

differentiating (7.4.4) once with respect to -f . 

Substituting \ /jL in (7,8.1) and (7.8.2) 

we obtain 

w I =° 

>=-« h-»fe 

that is, the case of straight annulus. 

7.9 HUMIHICAL EXAMPLE 

So far we have found the various flow parameters 
for arbitrary value of the ratio cr # jfow we proceed 

to discuss a particular case for which the radii-ratio is Vj., 
that is, 


cr ^ Va- 

so that 


(7.9.1) 


Lcr * * 0.6931 


(J^ a .) V * .480388 (7.9*2) 

3 

A * *1.082096, A » 1.17093, A * - 1.2670S7 

L, 

Lx. 

^3 


*0.608647 * , 
*3.700634 £ ) 
0.114647 P> , 
n.najsaa. & . 


t * -0.989286 X , 

^ * 0.048879 X 

f 

(a ■ 0. 116781 R , 

w • .000266 K . 


(7.9.3) 



f 


153 


H, 

* 1.020624 

0.028614 

R" 

Mi 

* 0.208266 

0.006700 


X 

* 0.058437 R' - 

2.166048 


7 

* 0.062600 

0.000605 

R. 


(7.9.4) 


Using these numerical values we obtain 

- 0.048367 * 3 

♦ 0.022644 Ai 

(7 .9 *S) 

- 0.097106 -f 2 " 

♦ 0.112718 -AM 

♦ 

- 0.219649 " (7.9.6) 

* (0.0686 * 0.000606 ft) ^ 

+(0.020302 R x - 0.812600)1 + (0.76-0.026864 R^) I* 

+ 0.007126 0.000046 ft + 0.000087 

- 0.000939 P-^ 7 ^) - 0.00804 

+ 0.161899 R" I 3 (^1) + (0.017834 0.641048) -s Jk-c 

- 0.000144 4 • + 0.006099 tO'-s^As) 1 ’ 

- 0.009899 a"-. 1 3 (A-S J r 0.006612 -5 (7*9.7 ) 


/ 

* 0.000266 ^ + 0.046671 -? 

+ 0.006804 - 0.003472 -? 

+ 0.047696 + 0.02444 L~{ 

- 0.073183 7* (AlJ 2 - 

t _ I 

Jp- ^ - 0.000266 ^ + 0.070110 

£ £ 
+ 0.061566 1 - 0.024306 -3 

- 0.003278 « + 0.024434 ^ 



_]L * -0.048879 ^ - 0.742968 -f + 1.635983 ts 3 

-0.833333 ♦ 1. 523143 -f 5 (4-jj - 0.49663 -f'C-k-O 

-0.585465 f(ivi) + -f (i-3+1.082095 luij ^ (7.9.8) 

7.10 PHDJkCflOK OF fHI SfRhAK L1H&S OK A HOiftAL dlCflOii 

Oeaa (1927} has drawn tine curves of intersection of 
the surfaces = constant with a normal section & = constant 
for the flow of an ordinary viscous liquid through a 
curved pipe and Jones (i960) has done this for the case of 
non-H ewtoni an liquids. We now do this for the present case* 
The curves have the polar equation 



(7.10.1) 


where R (which is independent of R ) is given by 
(7.4.3), and K is an arbitrary constant. Since k takes 
arbitrary numerical values, the relation between ~i and $ 
expressed by (7.10.1) depends neither on ^ n or on £ • 

For the case of CT -.A , (7.10.1) takes the fora 

= a- [o. 000266 * 0.045671 ~i - 0.048267 <5 3 

♦ 0.006804 * - 0.003472 -f + 0*022544 ) 


4- 0.047696 ■i (ta) ♦ 0*08444 

3 7- —I 

- 0.073183 ** ( br$) J (7.10.3) 


We shall use <J* ^ich is given by 
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^ ® 100000 cr' (7,10*3) 

Calculating the values of ‘V# (defined fey (7.9.6) ) 
for various values of f in the closed interval 0,6 $ f^i.O 
found that the function ^-//j attains its suliui value 
nearly at =» 0,80 (independently of Reynolds number). 

This means that ( a /r ) ( ) = 0 at -i * 0.80. Thus 

at the points ("fa 0.080, <^->rrr ), the velocity component 
V vanishes. For every particular value of a- , there is 
therefore a limiting surface 4- > (the constant 
depends on cr ) which takes fee degenerate form of a single 
circular streamline in a plane parallel to the central plane. 
The intersection of the stream line 4- - ^ with a normal 
section <9 * constant is a point. The stream line V-= 
is defined by cr' _ |0 /& when <r and the corresponding 

point of intersection with the section & * constant is shown 
in fig. (7.2) crV tP » for all values of cr , corresponds 
to the annulus-walls. For fee case of <r •=, Va_ y the curves 
of intersection of the surfaces corresponding to a-' *:3,3 t 4 f 6 f 10 
and with a normal section 6 * constant are indicated 
in fig. (7.2). 

7.11 VLOfeillt in m C^raAL flami 

The differential equation for fee stream line in fee 

central plane is 

__ , (>= ±s] (7.U.1) 

U w w 



which can be written as 


or 


ot<S> 

<=tA 


Vj 

[\-\h 


c = 




J 


c<c& 

cfi. 


. 'f 

A-c?- 4^ / 




+ M 5 . > 0 '. 

uhU 



l*' i?) 

(7.11.8) 


Equation (7.11*2) shows that the stream line in the central 
plane will depend on the ratio Si . Following Dean (1028) 
to a sufficient approximation, we may write (7.11*2) as 


-V <=(® _ 'i ^ A4»x-i ") 

'M ~ ' ' (7.11,3) 

the sign ;± being the sign of 4 , and the value of 4- -the 
one given by (7.4*4). Equation (7,11.3) shows that the flow 
lines will vary with a- and R . For cr W we have to 
evaluate the integral 


zie =. 


SACKS’ Lt ) . 

- ■ - ^ (7*11*4) 




where V - 7 is given by (7.9.5). 

The (a,^j relation given by (7.11.4) is seen to be 
independent of the ratio °-f jl • Ibr a given value of 
varles Inversely as R * The magnitude of &■ increases 







continuously with C and tend# to infinity as £ toads to 
unity* gineriea illustrations ar* s dw (lv«a for R » 100 
f" - \ *®* o* - 7 CXt Is doubtful whether the 

approximations made ia tbs above discussions are strictly 
valid for this oboist of f ), fits Integra ( 7 . 11 . 4 ) eta 
b« evauated numerically to giro 6 as a function of C * 

She results of such a numerical integration are given in 
table 7*3* 

fable 7.3 i Values of & in degrees for corresponding 
values of £ 

«ria^v^: ; KiAiiwsv'ar'io’t:.;::. *»Jc,va,vK-:,ik^f ::nA' .':sna^. : i> er.a.vCi'eeiusis aesorMa:'. m.a:vt. a-':mv . :-Kir:"-::ax ■; eetwi s: Mr*«e :r. ; av.;c::>:ia_. 1 i„:::. rtsL'rvsui..* eiuimnr iM ev saur * 

Values of & in degrees for corresponding value of C 

vaw:iic*owir - % 1 ^v»Ker:r*vir 

C 0*40 0*66 0.70 0.75 0*80 0*85 0.90 0*95 

<9 8*94 14*50 88.01 88.54 33*43 38.45 44*14 54.81 

qiseewiawjjuratBOuVK^^ ease. wirm. .%««**'. **ir »ww;,;v*iev:Ufta«s.r;.vc:;.i:^:;:KJ»-»sr^iy:/:;;ia',::rvi^aa:,wi'r«rui» 

fhe form of the stream line for the particular value 
cr- \ Is shown in ths figure (7.4). 

7.18 the irrsof of lAi&us uno oi raocxti imo 

In order to examine the effect of the radii ratio 
cr- ©a the velocity field§ we have performed the mmerica 
computation on equation* (7.5.8) and (7*3.10) - the pertur- 
bation velocity end the primary velocity lee the 

straight annulus. The results have been given in tabular 
fern, fables (7.5) to (7.10) glee tb# wworieal vanes cf 
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( 


_w_ 

A *' 1 


)• < 


) Old ( Kr 


) at various 


points < on the horizontal diameter of the annulus for fixed 
value of the Reynold's number R = 200 and for the curvature 


q /jl =-o 2 • the choice of R is arbitrary since it is within 
the limits of permissible Reynolds numbers since it is less 
than 2000 the critical nimber for stream line flow mentioned 
by Dean ( 1928) and also since ‘o- the secondary velocity is 
everywhere less than 5 percent of the primary velocity i • 
This is Justified for the first approximation taken for the 
problem* The radius ratio of the cylinders has been taken 
to vary from .06 to *1 at an interval of *0l . The contribu- 
tion of w on the convex side O ^ and the concave side 
- 9v) °f the horizontal diameter is Just the reverse of 
each other i.e. on the concave side up to the point ~f * *20, 
the contribution of this velocity is positive and after that 
it becomes negative. On the convex side it has Just the 

i 

opposite effect. The point where the perturbation velocity w 
changes sign move towards the outer wall as we increase the 
value of R and it reaches somewhere near S, * .36 for 

<r * .1 . 

In tables (7.11) to '7.13) we have tabulated the 
velocity fields for fixed value of the radii ratio <r *.026 
and a /> * #06 for different values of the Reynold number R 

at all points on the horizontal diameter. fbr R * 60, the 
contribution of w is positive at all points of the horisontal 



diameter on the concave side whereas It is negative o© the 
convex side. It Is also more prominent towards the inner 
vail than towards the outer wall, is R Increases to R « 76 f 
we see that up to f * .40 the sign of w is positive and 
then it becomes negative and again at ~i * .80 and onwards 
it becomes positive* this shows that for some combination s 
of R and cr the secondary velocity may be aero, llama R 
increases further to R * 100 , the point of change of 
sign move towards the inner wall. 
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Table 7.33 i Yeleeity field for R * 100, — * .OS 
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CuitV*C£l¥E EEA'JJ- XK aU H fhSL 1M a CUHVsO PlPs OlAUttfLUut, FLOW 

8.1 itimMcma 

Consider & laminar and steady flow of a Newtonian 
fluid through a curved pipe whose cross-section is a simply 
connected circular region and the curve Joining the centres 
of cross-sections Is a circular arc. Let the wall tesaperature 
be constant with respect to time as well as direction, fluid 
properties (i.e* viscosity coefficient, thermal conductivity 
coefficient etc*) are constant, free convection effects are 
negligible (i.e. flow is a forced flow), heat radiation and 
heat source (or sink) distributions (such as t chemical or 
nuclear etc.) are absent and both velocity and fcenperature 
profiles are fully developed* We confine ourselves to the 




case of liquids only so that we can omit heat of compression. 
thU8» dissipation faction Is the only representative of 
heat transfer in the problem. The problem for straight duct 
was first of all studied by sehlichting (1951}. The solution 
of she corresponding problem of straight pipe of circular 
cross-section can be seen in Landau and Llfshits (1950). Ihe 
purpose of this chapter is to extend the work of Dean (1927) 
to Include the ease of heat generation due to viscous 
dissipation. We follow Dean and evaluate the temperature 
profile at the same level of dissipation at which the velocity 
profile has been evaluated by Dean. 

8.2 aii4iiii4L mn tti&n&nw&i bimlv 

A circular pipe of radius a is coiled in a circle 
of radius t » i.e* the locus of the centres of all cross- 
sac t ions makes a circle of radius • The coordinate 
system used by us is the same as was adopted by Dean (1927) 
and Jones (19b0)in discussing the curved pipe flow problem 
and is represented in the figure (8.1). O V is the axis of 
the circle in which the pipe is coiled. C is the centre 
of the circular cross-sect iun by a pi a* which passes through 
Oi? and makes an angle 8 with a fixed plane o) passing 
through O'i? * C 0 is perpendicular to O v and of length £ • 
The plane through O perpendicular to O V will be called 
the 1 central plane* of the pipe md the circle traced out 
by C its 'central line*. Any point P of the cross-section 


is referred to by orthogonal coordinates 'Yj 4 s > <9 where V 
is tbs distance C- P sad ^ is the angle made by c P with 
th« lias through C parallel to OV , fhs surfscs of tbs 
pip# is given by Y— cu , a being the radius of say section. 
Th# components of velocity corresponding to these coordinate# 
are C_ u > VY W)j, U i§ therefore in the direction C. P_, V 
perpendicular to U and in the plus# of the cross-section, 
and W perpendicular to this plane. The general direction 
of flow will hi taken to be in the direction in Which 9 
Increases. 

When the incompressible Newtonian fluid characterised 
by the rheological equation of state 


tcj - ~~ ^ 


(8.2. l) 


flows through the curved pipe (due to the fall in pressure 
along the pipe) with a uniform will temperature, the tempera- 
ture profile is changed fro* the uniform profile at the 
entrance to the characteristic one down stream and at a 
later instant a fully developed temperature profile is 
formed inside the pipe. To this temperature is added the 
heat of viscous dissipation. The viscous dissipation function 
for the above liquid is given by 


4 = Zfx. ey e>j , -N<«.s.*> 

, 2 + «? +* t&r 2 v* 2 $J,J 



m 


the general energy aquation, for the liquid# of constant 
properties with zero internal heat source and sere radiation 
it 



where K is the t hemal conductivity, f 3 is the density end 
Cj, is the specific heat at eon st ant pressure, the quantity 



is usually very small for liquids. 


It is, there- 


fore, taken equal to zero in the present case, fhe equation 
of continuity is 


\j i j __ j_ ~&v 


following Dean (192?) we assume 


U - U, 


V- ^ 


_ A~ ^ 

~ 4/0. 


and 


-4-B + i 3 


where 14,, ^ and I 3 are referred to the secondary 

flo„ and all ar. tak«n to b. -all of to. ord„ of f ^ 
u u U5- are independent of & « Ms further rewrite 


(a>. 


b 

1 I 


/L. r^t) 

4m- v y 

-/l*. + 6 - 


(8.3.6) 
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wbiek are nothing but the velocity and the pressure fields 
for the flow in a straight circular pipe under the constant 
pressure gradient. Thus in this case we can regard ta , UJL - 
and Ip as snail. With this the equation of continuity 
takes the fora 

^ t + L. It _ 0 

-v Y ^ (8.2.6) 

and this suggests the form for the stream function ^ as 

LU^ _j _1± and i*=|£ . (8.2.7) 

~i ~2i hr 

Thus the strain-rate and the stress components are simplified 
and the simplified equations of motion reduce to 




2 ± 

Or 


At 2 r 3 

V V< : £i_ r Sy2- 


L?* _u.it 1 

+ y 3r 4 'r^4> 3 -j 


(8.2.8) 
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T. 

3 V 7 


Jj(®< 


2 . 9 ) 


p A_^t + M = A-cfe <«•*•*» 

r ^/U. h<p %<- 


where sms found to bo constant and has been taken to 

^ 2 


be sere for simplicity 



8*3 MN£fl&¥ S^IAHOMS 

&ot the straight circular pipe, we hare 

(8.3.1) 

o = ^ v. - A r - + ^7 v J (8.3.8) 

»fi<4 tr £ = f- _.^ Z . ( cb-^) (8.3.3) 

where la the constant wall temperature, and ^hl is 

the temperature at a point inside the straight pipe. Just 
as in the case of velocity field in the curved pipe (see 
Dean (1987)), ve write down the teaperature field as a 
perturbation over the straight pipe case teaperature 
and this perturbation is small due to the large radius of 
curvature L % let 

TiC*,*) = t/eCK*) + 

or siaply (8.3.4) 

r e = t IL +■ t L 

Here occuring in (8.3.4) has been written at the level 

of approximation at whieh the velocity components and pressure 
have been written in Seas (1907). 


fbr liquids, the energy equation (8.2.3) is written 


AS 


f c b U- 4 tjt 
r Ay 
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Ji Y 


USi^^P \ 

' "L ) 
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£ civ 


+ k (_ f -L 3 ^ + 

^ "Y ' Y Y 


j ^ \ 

v L ^“4^ j 


(8.3.5) 


Sine# 


cLbtt 
cL y 


yy 

/6M~k 


and 


<£~r 


Ay_ 

2 M- 


(8.3.6) 


Equation (8.3. S) simplifies to 



/6/Y /c. 




-f 


zl 

Y 2Y 


4 -^—- ct-S~ ~r ) 
/6/LLl^- y ' 



i^l + j. §^1 j_ 1 tt ) 

2t Y 2 - Y 2 y + T 1 - / 


(8.3.7) 


8.4 aoiunoiid 

following Dean (1827) and Jones (18oO), the equations 
of motion (8.2.8) to (8.2.10) can be solved by assuming *fO ^ 
and V of tbs fora 


s*id 


i b .= ^ 4> 

r bp- L 


Uy 


l 

p ^ u 2 c 4>j 

~HP-i 


\ 



bti) 


(8.4.1) 



where j and \> ) are dimensionless functions 

X -y*» 

of ^ " X" * ^ bUs equations* when substituted in the 

equations of notion (8.2.8) to (8«3 a 10) reduce to 


-kt<- 4 2- d - 


(8*4*2) 

- AO _4/ = 

-U <d- A lb' 

* & r 

(8,4.3) 

— 3. P -f & — 

A w' 

(8.4.4) 

vtMr* j 

R - P Aa - 

p ^ 

] r *^= — 


and the operator is defined by 

. !L _ i- 

?V + -i 4 " 


Consistent with the 

fora of <** and 

f ^ t**e 

solution for will be of 

the fora 


* 2. 5“ 

f A a* 


(8.4.5) 


2 ¥ M 

Substituting the values of ^ end 

we get 


etc. in (8*3*8} 


k P y + /6 ^ Vt ~ A t(L 


(8.4.6) 



Ql»«n this equation n „an date nine t' e provided 
and are knovn and these m have froa Jones and Dean 


as 




ft 


£44 -<?4 3 +- ^4 i - 4' 
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(8.4.7) 


w'= -f£4-4 J + £l , 

4 U tfSao 1 - 15 V (8.d.8) 


and 


(='= 0-74 


(8.4.9) 


thus 


zl t 


(2- - 




- — S3L ^ ~ C) - ^ c^f 4 - 4 3 j 


izt£-*C) + (i^ - ao^\,sc 4 - 7*4? 4^J 


(8.4.10) 


and t her® fora 
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3^ 


(4 t t ) =_ 

^ ^ /44o 




2> a 


+ „A^_ £ (jqo -&(>%+!> oo£> -/7r4' Wl)+ 8,4 


(8.4. U) 


6> i heing arbitrary eons tent of integration, this equation 
again on integration giro* 



1st 
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i. _ £‘p T 4‘ 


mare 


/ 4 o - 4r4'V / 4 & j 


".IL 


+ -zoo&i*^ 4 ~s^4 6 *-s4 i? ) .+ _ 6 '^ + £2 

<8,4.13} 


; -f- ~C +- 


vh«r« i g again an arbitrary cent t ant of integration. If 

tbia constant 8^ were finite and hetce the temperature 
would be infinite at the centre of the cross"- sect ion. since 
the temperature at any point inside the pipe has to be 
finite and not infinite}, the constant &?_ must be zero, 
in order to determine i3 ; , we have the thermal boundary 
conditions t t - = o at 4 = ^ froa the prescribed condition 

that T g - fc-oj at 4-^ and therefore froa equations 
(8.3.3) and (8.3.4) ve get 


3, 


x + "*1* + Rl 

3 x !2co 


and therefore) we get 


(8.4.13) 


= C 7: + 




7 , //& fir t A 


6 /yzS Yo ^ Lid0 


gjPr 'if C) 

\lZ%o v - 


6 , 


4 - . 6 *" - . r'/qS'-ZLOoC+m £-1^ + 3>tf ) (8.4.34) 

%2%-ao L y 



&osi equation (8.4.14), we note that ^ =• ° when 
4 = ° * A* a natter of fact *4 should vanish at the 
centra of th« cross-section, otherwise the solution will not 
be of any physical significance, this is explained as follows* 
Consider two points on oc , one on the left hand and the 
other on the right hand sides of c . at the left hand 
side point 4> - and , while at the other J 

and i^,c£ = / * There will therefore be a finite difference 
in the values of (and hence, in the teaperatures also) 
at these polntsi however closer the points Bay be unless 

^£-0 at 4 - o * This iaplies that the temperature is 

/ 

undefined or discontinuous at 4 = o when it 4= c * 

This Is absurd from the point of view of physics. Thus, for 
the solutions to be of any physical significance, *= c *ust 
necessarily vanish at 4 = o * 

Thus the tcaperature field at my point inside 
the pipe is given by 
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Tbs rat® of beat* transfer 


? ic 

44 ' 


is givsa by 


4 ^ j O- j 


ft p r / - // 1- 4 - 3 / s- 4 6 +i 6 s-£? ~ 3 .z 4 ° ) 

/ 7 i’*a ^ y 


4 - yy %s£-idoo 4 + ^ 4 -3at 4 ■+ 33 4 3j (® ,4. ^ ) 


and Its value at tbs surface of tbs pips i.o. ^ - \ is given 
by 


I = 
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Figure 8-3- Thi ferfiperoture field in* id* the Curved pipe 
for Pr • 10 » ft *fS ah ail the radii vector*. 
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W# have earned oat numerical evaluation of the 
temperature field and toe local heat transfer at the wall of 
the pipe on I at 1680 Computer. The Prandtl number TV hat 
be®a taken to he 1 or 7.6 . The Beynold»a numbers R has 
bemi taken to be 60 or 76 or 100 along with the ratio of 
the radius of toe cross-section of toe pipe to that of the 
radius of curvature of the coil to be .01. tie have drawn a 
number of graphs showing the behaviour of the temperature 
curves and the local heat-transfer rate curves for various 
combinations of R and Rr » fire* the graphs we stay 
draw the following conclusions. 

(i) Graph (8.3), shows the teeperature field inside the 
pipe for Rr *1.0 and R * 60 on all the radii vectors. Me 
notice that the temperature inside the pipe on the convex 
half of the pipe is a little higher than toe corresponding 
points on the concave half of toe pipe. On the horisontal 
diameter on toe concave side at about 4 * .3 the teeperature 
attains maxima# value which is above the wall teeperature. 

(ii) Graph (8.3) has been drawa for R * 78 keeping Pr 
the same as 1.0 • Behaviour of the temperature curve remains 
the same but for a little change in temperature towards 
higher side. The point 4 * .2 also has a tendency to shift 
outwards* at no point m toe concave side, the teeperature 
inside pipe la above the teeperature of the wall of tho 
pipe* 
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(lii) Oripb (8*4) shows the temperature curve for ^ * lOO 
and remaining the saae at 1*0 • Here also the behaviour 

of the temperature field remains exhibiting a more appreci- 
able difference in temperature on various radii vectors* At 
corresponding points on the convex and the concave side ef 
the pipe, the temperature difference goes on increasing till 
4 * aod then goes on decreasing mid becomes the same 
*t 4 • 1.00 as la every other case. 

(iv) Graph (8,8) shows the temperature curves ^ * SO and 
Rr * 7*S to illustrate the effect ©f increase of ^ * 

Behaviour of these curves is again the same, mid there is 
appreciable increase in the temperature at a point for the 
corresponding curve when £ * SO and K- * 1*0 • On 
*= I “ d *---1"* Sit f.r.no. 1. maimm 

at about 4 * *0 • 

(v) Graph (8.6) shows the temperature curves for ^ • 75 

and & * 7.5 on horisontal diameter, the vertical diameter 
and the radii vectors <f> - f end <p - - f • The temperature 
on the vertical diameter remains the same as in all the other 
cases, end this goes on increasing with increase in R , such 
that on the horisontal diameter on the convex side of the 
pipe, it attains maximum value at about 4 « .4 * The maxi- 

mum temperature difference is observed to be on the points 

4 • *6 on horisontal diameter on opposite sides ©f the 

centre of the pipe* 
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F ioure 6?‘The temperature field inside the curved pipe 
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Figure. 8-8- The wail-beat transfer rate of the curved pipe at 
various points on the cross section of the pipe for 
Pr = 75,R=75. 
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(▼i) Graph (8*T) shows the teipintura curves for ^ * 100 
K * 7.5 on the same radii vectors as before. Here we 
notice that the temperature difference at about ^3 a ,46 
still increases with the in R- and R .On the horizontal 
diameter at about <4 a .§*» the temperature difference is 
maximum on two sides of the centre of the pipe. 

Since the expression for the skin friction does not 
involve R , we can say that the higher temperature differe- 
nce in the liquid on both sides of the vertical diameter i# 
more due to the value of R than due to the rise in eklm 
friction, because the skin friction is higher on the convex 
side of the pipe than on the Gcneave side. 

(vii) Finally graph no. (8.8) shows the rate of heat-transfer 
on the surface of the pipe at various points on the cro se- 
lection of the pipe* He observe that for all combinations 
of R and K- , the local heat-transfer rate is highest 
at the end of the horisontal diameter on the convex side and 
lowest on the other «nd of it. it the ends of the vertical 
diameter, it is constant and independent of the radius of R 
and R • the rate of heat-transfer goes on deereaslng as 
we go on deereasing either R or R . 

fhe effect of curvature ©a the heat transfer is 
larger t han that on the skin frietion m&d this is expected 
because the skin friction 1« proportional to the veleeity 
gradient end the local heat-trenefer rate is proportional 
to - the second power of the ooleeit|^jtred ion <»* 

' A 


CHAPTER II 


CONVECTIVE HEAT TRANSFER IN A CUBVES ANNULUd SfRlAWLIrfE PLOW 
9*1 INTRODUCTION 

’ In the seventh chapter we extended the work of Dean 
<1927) to include the case of curved annulus geometry and 
confined ourselves to the analysis of hydrodynamic problem 
only. In the preceeding chapter, we extended Dean’s problem 
to Include the case of convective heat-transfer due to viscous 
dissipation in curved pipe flow. In the present chapter we 
wish to extend the study of chapter VIII to include the case 
of curved annulus geometry already considered in chapter HI* 
In other words the purpose of this chapter is to study the 
convective heat transfer due to viscous dissipation in the 
fluid to t ion considered in chapter VII. Wore specifically* 



our purpose is to consider the teat generation due to viscosity 
in laminar, steady, constant property fluid and fully developed 
(both hydrodynamic ally and thermally) forced convection of a 
Newtonian liquid inside a curved annulus* the temperatures 
of both the inner and the outer cylinderical surfaces are 
constant (at the same temperature) with respect to time and 
space coordinates* This is only for simplicity that we have 
taken the temperatures of both the cylindrical surfaces to be 
equal* The problem of curved annulus heat transfer is of 
interest in double pipe heat exchangers and in nuclear 
reactors design* 

9.2 MATHEMATICAL ANALYSIS 

Here also we adopt the coordinate system given in 
figure (7*1) in chapter VII. The assumption that Q/jg is 
of first order of smallness (i.e. (<*/*) 2 5 (« /4 j i . -- • are 

negligible as in chapter VII) is also retained for the 
purpose of satisfying energy equation in the following 
analysis* 

Mte know that the velocity profile for a fully 
developed forced flow under a constant pressure gradient 
of a Newtonian fluid along the axis of the straight annulus 
is given by 

w ' - A*r ( (9.3.1) 

where „ A = - , " cr ~ , h and A is the constant 



pressure gradient along the axis of the annulus , and 


^ - — A-t_ B 


(9i 2*8) 


together with Vj^o and v* o. 

The energy equation in the ease of straight pipe cited 
in the previous chapter holds good for straight annulus also* 
Solving that energy equation under the condition tiz ^ ^ 
at and = cr , we find that the teaperature at any 

point inside the straight annulus is given ty 


'IJ2. 


-t- AV* r l+^A+^])Kt -s.A % ] (9.2.3) 


L 


wh.re t = -tv—' ’)] 

and U is the teaperature at both the walls. 

Since we have assumed the curvature of the annulus 
to be very small, we can suppose 

T = t u -+ t* (9.2.6) 

similar to what we did for velocity components and pressure 
in chapter VII. T , given by this equation has been written 
at the same level of approximation at which velocity components 
a n d pressure has been in the preeeeding chapter* Bmp# 
is the contribution to the teaperature field due to the 



curvature of the annulus, and hence forth we will call this 
as perturbation temperature* 

ihe physical picture of the energy equation suggest 
that the perturbation temperature t*. (on the same lines as 
velocity, stream function and pressure) should have the 
form 


t t 


. .A 1 ^. , t H) . + 


(9*2*6) 


along with the forms of to, 4- and jb as assumed in tide 
VII chapter* 

Equation (9.2*6) gives 


32 

2 4 L 


= _ tl l-S ) . £uXq> 


(9.2*7) 


Thus the energy equation ( 8 .^ 3 ) becomes 


P c l> ^ 


cL 61X 

a-i' 


( 

V 'SsC 


_]_ ^JL , J_ . \ 

■S St s x B-4* ) 


j- Cui 4_ <=*£,£ - K <rfW| ( ctw _ t^! £44. \ 

4 c<-r r V a-i Q ~ J 


The solution of the above equation will give the 
temperature , knowing (9*2*3) *“ or ***• straight 

annulus* 
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9.3 SOLUTION FOB £* 

Substituting the values of to , and w the velocities 
for the straight annulus and its perturbation velocity 
respectively for the curved annulus from chapter ¥11 we 
reduce (9* 2.8) to the form 


p .( • ^'1 

" I” M s 


= k ( + -L JL 6- ^ h . ^4 

i ”*1 “-P- e / x.” ' 


(9.3.1) 


which further reduces to 
( . X i l' \ _ 

v A-ii ^ -t XA ~ c ' " 


-a p,. r ^ ) c/'^r- +#>4 

+ s(a-i+Aj^' -v [_ -i 

_ AD + ^ X + lA + 1 1 A j 


A 


where the non-dimensional quantities are defined by 

3 - *3T 

Pv - 

prandtl number - “ k 


(9.3.2) 


pAs 


Beynold's number R = 

AV 

Thermal conductivity *■ = 

tbm «quatH» (9.3.8) itl.ll «i»lT»a X<5 aad w'C-« ) , 


(9#3*3) 



and these we have already obtained in chapter HI and are 
given by (7,4.3) and (7.5.2) respectively. 

Substituting all these values froa the VII chapter 
i# (9*3*2) &s 




” o •+“ -v- *T -v* *f 

3 * 1 


7 <J 3 

+ C.fe ~S ■+ c 7 ~s ■+■ c 1 5 -"S ■ ) 



+ c,* v+ C„s ) M 


) 


(9.3,4) 


where c, 3 t,g are given by 


C n = 

^3 = 

- 
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C'C, = 

c 7 = 

C.& = 

Cc, = 
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_ f . 2- 
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R-VS 

— 

~~ S ‘ q * -+ ut> ^ 



= 

— S/t^Qg -v sa bm 


C-IS 

=: 

S (> q 7 - 3) <<-cj j 4~ & ( 3 A b 1 3 -*- X b 1 H ^ 



= 

*i_ 

— s ( A 9 C -4- 2.A O.^ ) + 8 ( c b, s 4-S-A b, a ) 


C »7 

=■ 

- S + So b lx 


C-i* 

=r 

i. 

- sa a s 

(9.3.5> 

where s = 

. A?*- and the constants a l) c il <3^ and 

ki ,K -- 

ar® given by 





h /R *, = — V<t 8 
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0<3- F ) AR • ^7 - C^A +S-J^’ 1 )/3i 


q 3 

= 

(--p - 2-<A -I7> 1 )//a8 j qj - f AR 


Ql, 
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(9.3,7) 


Xhe final integration of equation (9.3.4) yields 


tjlC-S) = + ( i±' -4-d 1 )j + d^t 3 + ^ 3> 'f +d h -f +d s S 

4 . (^^2. 4 . cfg-j 4 - cAq 4- =Jlo -£ -V ctw-i -+c<U~f J iMf 

3 £" ^ x. 
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(9.3.8) 

wliere and *V are constants of integration and ct > > j 

are given by 

d i = (-SCj.4 e S- c t 4 )/<S’ 
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4n = 
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4\^ = 
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dm = 



(9.3.9) 


fhe constants of integration N» and Na. are to be determined 
from the boundary conditions 


£*( i) 


= o 


and 


£jl(v) " O 


(9.3*10) 


thus we get 
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j (9,3* 3*> 


Ttius knowing the perturbation temperature with 

the help of equations (9.2,6), (9.3.8), (9.3.11) and ( 9 . 3 . 1 g) t 
we get the temperature at any point on a cross** section of 


the curved annulus by a plane & * constant for fully devele— ■ 
ped temperature profile case, as follows * 


"T = t u + (* 


too + r 1+ ^ m h ln-j - T- j 

6^ mk l 


+jAL. tlc-S) .^4 


CMP K1 


which may further be written as 


T - too 


[ 1 + + - bx-s^-s* ] 


+ A . L(-0 ,jStU4 


(9.3.14) 


The second ter® on the right hand side of the above eqmet&mm 
(9.3.14) is to be regarded as the super-imposed temperature 
over the temperature in a straight annulus ( the first tesn 
on the R.H.8. of (9.3.14)) due to the curvature in the annulus. 


This will be zero for large radius of curvature i l.e, for 
the straight annulus* 

9*4 IKE HEAT-tRASSFEB-fUrE Af THE WALLS OF THE AHNQLUS 

To find the heat-transfer-rate , we differentiate 
(9*3*8) and get 

^ °l I + c ( 8 ) ( i ^cl'X.-VcAc i )-J 2 ' 4 .(S-c(3+.dlo)-f tt 

et-S S* ' ~ 

4 & lo 

+ (7^ +.oln) i -v -v d\r)^ -\r u d fe s 

+ j (_ a<=iii _ a 7 j/-s^ + (_ J s •+■ csotq -*-x<=k\ S 

i ^ O ® 1 

■+■ (^Sc^\ 0 -t 2.c* ic.')'S -v(7<=A||-+-i Q <\7^-i-+ c 1‘H.ij.-5 ^ ( Ln.'f J 

r <=*-V 3 . V 

• Jr l - -v ) -v- ( ScA,^- +3<=<i^ ) -i 

H C? , *- 

+ S"o<rg 1 1 -t- V <=< (7 tT j (lu-j) 

{ °*> 8 -V Scl,^ -i 2 ' ^ (9*4,1) 

This is regarded as the perturbation ter® fro® the heat- 
transfer rate over the straight annulus and ue get the 
following expression for the heat-transfer-rate at the walls 
-<r and ~i - 1 (as the case may be) for the curved 


annulus 


cXT 

& \ kk _ _ 8 A i 

C^aV^Mk) L “S J S=cr 5 d 

-v- 3_ • f^ a ' |.^4 

'iso-j ^ (9#4#8) 


9*5 DISCUSSIONS 

We have performed numerical computations on equation 
(9*3*14) on IBH 7044 Computer* The numerical values of the 
temperature inside the annular region are given in tables 
(9*1) to (9*5). Specifically these tables given the values 
of (T- boo) /6<<hk) at several points on both sides of 
the horizontal diameter* As in chapter VII the choice of 
^ - the Reynold *s number Is arbitrary in the range of 100 
to 500* The radii ratio of the two cylinders has been 
taken to vary from ,1 to ,6 . The Prandtl number Pv has 
been taken to be 1, and it is also in the range of the 
Mewtonian liquids. The curvature of the coll has been taken 
to be *03 which is sufficiently small in comparison with the 
radii ratios. Apart from the points on the horizontal 
diameter we have also calculated the temperature profile on 
the radii vectors = o and <3> = ± ^ . In general we 
see that she contribution of the perturbation temperature 
tji on the concave and eonveat side of the vertical diameter 
is just opposite in sign. Thus we discuss all the cases 
on only the convex side <3> = 3E of the horizontal diameter* 





Table (9*1) gives the temperature profile t|.g. for 
the straight annulus and the perturbation temperature (■£ over 
assumed due to the curvature of the coll* The 
Reynold's number has been taken at 100 and the radius ratio 
« 1 the curvature at *02. We see that in case of straight 
annulus, the temperature rises more steeply near the inner 
wall than near the outer wall* In between the walls it 
becomes more stationary with the maximum at t * *46. Hie 
also find that the Reynold number R has no effect on this 
temperature profile for the given radius ratio of the two 
concentric pipes* Hie perturbation temperature increases 
in the beginning to "f * *25 and than it becomes negative 
till t * *40 and then again it becomes positive and remains 
so till 'f * .76 to become negative again till the outer wall* 
The net result being the resultant of the two to make the 
temperature more stationary in between the walls* Other 
things remaining the seme if w» increase Reynolds number to 
200, the temperature profile for the straight annulus remains 
the same but the perburbation temperature begun s with a 
positive value from near the inner wall and remains so till 
* .70 and after that it becomes negative and remains so 
till the outer wall* This results in a steep fall near the 
outer wall. As we increase R the contribution of the 
perturbation temperature becomes more mad more negative 
after the same point ~i * *76. This tendency of the tempera- 
ture profile remains at the radii vector ^ also. 


4 U 


be coming zero at the vertical diameter* Thus we can uy that 
the effect of is maximum at the horizontal disaster and 

becomes smaller and dialler as ve aove upwards along on arc 
somewhat parallel to the outer circumference reducing to a 
zero at the vertical diameter* This happens for all R 
increasing* 

fable (9*2) gives the temperature profiles in the 
same case as above but the radii ratio increases to <r * .g. 
We see that for R « 100, the perturbation temperature is 
everywhere negative on the horizontal dWdsxoa the convex 
side of it* When R increases to 200, all other parameters 
remaining the same, the temperature fc* becomes positive 
in the neighbourhood of the inner wall (-1 * .26 ) and 
Immediately afterwards it becomes negative and remains so 
throughout upto the outer wall* For Increasing R the 
tendency of the point of change of sign is to move towards 
the outer wall reaching ~t * *30 for R * SOG. For value# 
of cr ^ and R ^ 3oo the perturbation temperature 
remains negative throughout the convex side of the annulus, 
with the result that the temperature on this side is below 
the wall temperature whereas on the concave side it is 
higher than the wall temperature* This is evident from 
the tables (9.3) to (9*6). Apart from tbs tables ve have 
also found that no change of sign for the perturbation 
temperature occures for lower Reynolds number is the range 
of 26 to 76. 
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Jbi some values of stir swcrO) tlia gpi rtuftatiou 
temperature values are negative at sow paimnisi, and positive 
at other points on the horizontal dl.a©«ts«r, f For larger 
values of cr the values of the sms© tmspsrsatwre are 
negative at all points of the saute <3 lessee ta? , this further 
shows that for large values ot O' tlae tteei^rinture la 
decreased on the convex side and lncr«eaes»«] outhe concave 
side* 

ye have also computed the local! Wan aU«Ursnsfer-rate 

from equation (7.4,2) on the walls cf tlfae muatilus for 

- “ #■. 

the same set of parameters as before, **itn« poointa on the 
surfaces have teen taken to be the ra_dizi yoettors 4> = ° , 
±]£ ^ ± £ . Here we set that at fc-h emdJaof the vertical 
diameter i.e. 4> =• o the pertur hatLoin t.empOBMture is always 
zero for all set of values of P* am* cr nlltfrcause of the 
multiplication of ^4 . I he perfcur-bsatriCT temperatures 

at <4 = ± have the same nnaerlcaal valwa hut have 
opposite signs, the same is the mane at *4 - ± ^ hence 
we study the perturbation tempoxatur© onlly i-J the ends 
of the radii vectors 4* * -t+ suk* *4— - 

Xhe wall- he at -transfer-cate straight annulus 

is not affected by the Reynolds nm!m: ^ »cr ft* prandtl 
number Pv t but the wall«heat“traa_sf«er*Tritci is of the 
curved annulus is affected by ’them, ®ie* tsjfcdbl* number 
gives the heat-transfer-rate at tx> the ta© »naar and the outer 
walls of the curved annulus for fined imILmm of the radii 



ratio ^ 1 % the Praadtl number Pv=sl ‘0 sajj the curvature 

q 

1 i,o!L bat for different values of R . from this 
table v® find that as R increases tbs stjoitudo of tbs 
secondary beat transfer rats at both the nails increases 
whereas the signs of the heat-transfer rate at the inner wall 
being negative but that on the outer mil is positive* In 
other words we can say that due to the curvature in tide 
annulus the heat on the inner wall flows towards the wall 
and that of the outer wall moves towards the liquid, fhe 
net result is that upto R * 300 the total heat transfer 
rate on the outer wall of the annulus remains negative and 
after this it becomes positive there also. 

The effect of cr on the wall-heat- transfer rate 
has been tabulated in table no* 9.7 . Here we have taken 
the same parameters as before for fixed Reynolds number 
R ** 100 but for changing cr • We find that as cr increases 
the secondary heat-transfer-rate at both the walls increases 
in magnitude, for cr ^ « 1 » the secondary-heat-transfer- 
rates at both the walls arecpposite in sign, but for larger 
values of cr the two rates are of the same sign, that is 
to say it is positive on the convex side and negative on 
the concave side for all values of 4> • We also find that 

the sigh of the second ary-heat-transfer-rate changes for 
some values of cr in the vicinity of cr= -2. * the net 
result in this case Is that the total heat— transfer— rate 



on the inner wall for all values of 4> remains positive 
on the convex side whereas on the outer wall it changes from 
negative to positive as cr . increases* 
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9*3 i Tfa# t«ap»rafcur« profiles on 4 tor ^ • 1*6, 
1 m 306, j- m .02 aad cr * .1 
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